o 
o 

(N 
< 

Angle and Volume Studies in Quantized Space 

> 

o 

oo 

Q ■ Swarthmore College 



X 



Michael Seifert 



March 19, 2000 



Abstract 



The search for a quantum theory of gravity is one of the major chaUenges facing 
theoretical physics today. While no complete theory exists, a promising avenue 
of research is the loop quantum gravity approach. In this approach, quantum 
states arc represented by spin networks, essentially graphs with weighted edges. 
Since general relativity predicts the structure of space, any quantum theory of 
gravity must do so as well; thus, "spatial observables" such as area, volume, and 
angle are given by the eigenvalues of Hcrmitian operators on the spin network 
states. We present results obtained in our investigations of the angle and volume 
operators, two operators which act on the vertices of spin networks. We find 
that the minimum observable angle is inversely proportional to the square root 
of the total spin of the vertex, a fairly slow decrease to zero. We also present 
numerical results indicating that the angle operator can reproduce the classical 
angle distribution. The volume operator is significantly harder to investigate 
analytically; however, we present analytical and numerical results indicating 
that the volume of a region scales as the 3/2 power of its bounding surface, 
which corresponds to the classical model of space. 
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1 Introduction 



Most scientists would probably consider the two "great theories" of twentieth- 
century physics to be Einstein's theory of general relativity and the theory of 
quantum mechanics pioneered by Bohr, Heisenberg, and Schrodinger. Gen- 
eral relativity describes the interaction of matter and energy with the fabric of 
spacetime (we experience this interaction as gravity), and accurately predicts 
the large-scale structure of the Universe. Quantum mechanics, on the other 
hand, deals with small-scale interactions between matter and energy, success- 
fully explaining the properties of electron, atoms, and other minutia. 

General relativity, in some sense, could be considered as the last of the 
classical theories of physics: it considers space and time to be continuous, and, 
with sufficient information about the current state of a system, allows us to 
predict the results of any future experiments with certainty. This is at odds 
with quantum mechanics, which often requires observable quantities to have 
discrete spectra; which says (by the Uncertainty Principle) that one's knowledge 
of complementary variables will always have a certain degree of uncertainty — 
we simply cannot have complete information about the current state of a system; 
and which only allows us to predict the probabilities of experimental results. 

One would hope, then, that there would be some way of reconciling these 
two world- views. The program pursued by theoretical physicists since the estab- 
lishment of quantum mechanics has been, given a classical theory, to attempt to 
construct a quantum theory which has the known classical theory as its large- 
scale limit; this process is commonly known as quantization. In the case of 
electrodynamics, the other major classical field theory, this program has been 
successful: the theory of quantum electrodynamics (also known as QED) was 
formulated in the 1950s by Feynman, Schwingcr, and Tomonaga. Out of QED 
grew the Standard Model, which includes not only electromagnetism but also the 
strong and weak nuclear forces (which have no corresponding classical theories). 

However, there does not yet exist a consistent quantum theory of gravity. 
This is partly due to a dearth of experimental data around which to build 
a theory. To obtain data to help us in this quest, we would have to devise 
an experiment which would (directly or indirectly) measure the gravitational 
interaction on a small scale. This pursuit is stymied, however, by the incredible 
weakness of the gravitational force compared to the other three forces of the 
Standard Model on those scales. An illustration of the vast difference between 
the strength of the electrostatic force and the gravitational force is attributed to 
Feynman: If you jumped off the Empire State Building, the gravitational force 
between you and the Earth would accc;lcratc you to a speed of about 80 m/s over 
the 350-metre height of the building. However, when you hit the ground, the 
electrostatic repulsion between the molecules in your body and the molecules in 
the sidewalk would decelerate you an equal amount over less than a centimetre 
— a minuscule distance. 

Moreover, the scales on which quantum gravitational effects might possibly 
be observed are minuscule. The fundamental constants governing a quantum 
theory of gravity are h (Planck's constant), c (the speed of light), and G (the 
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gravitational constant.) If we express the units of each of these constants in 
terms of length, mass, and time, we obtain: 



(mass) (length)^ 



time 



length 



c 



time 



G 



(length)^ 



(mass) (time) ^ 



Playing around with these constants, we discover that the only way to create a 
quantity with dimensions of length is 



This length, known as the Planck length, is the scale on which most researchers 
expect to observe quantum gravitational effects. This is an incredibly minuscule 
distance; as a basis for comparison, the radius of a proton is on the order of 
10~^^ metres, a full twenty orders of magnitude larger! This fact adds to the 
difficulties in finding any experimental data to guide the search for a quantum 
theory of gravity. 

Undeterred by this lack of physical data, many physicists have proposed 
theories which attempt to reconcile quantum mechanics with general relativity. 
Some of these theories (which will not concern us here) attempt to expand or 
modify the Standard Model to include gravity; these approaches include the 
popular theories collectively known as string theory. Other attempts to recon- 
cile quantum mechanics and general relativity are more conservative, merely 
attempting to construct a theory of quantum gravity without worrying about 
the forces included in the Standard Model (at least for the moment.) The 
theory we will be dealing with in this work falls into the latter category: it 
attempts only to create a quantum theory with general relativity as its large- 
scale ("classical") limit. This is the loop quantum gravity approach, in which 
the fundamental variables of general relativity are encoded in loops (essentially 
closed paths through space) ; the observables of this theory then depend only on 
these loops. 

Since general relativity makes predictions about the structure of space, we 
would expect that this theory would make predictions about the structure of 
space on a small scale.^ In particular, this theory makes predictions about the 
possible values one could observe for such "spatial observables" as areas, angles, 

^ The reader may note here that general relativity does not just predict the structure 
of space, but the structure of spacetime. However, the theory of loop quantum gravity (as 
we will consider it) is a canonical theory, i.e. it requires a separation of the kinematics and 
the dynamics of the system. As a result, we will only consider the states of the system 
(corresponding to space only), and not their evolution. We must also point out here that while 
some theories of the dynamics of spin networks have been put forward, their time-evolution 
is still not fully understood. 




(1.1) 
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lengths, and volumes. In this work, we will describe recent research into to the 
spectra of operators corresponding to angle and volume. 

Note that this work is intended to be understandable to the reader without 
any significant background in general relativity. As a result, a large portion of 
the initial theory (Section ^) is couched in terms of a theory the average reader 
may be more familiar with, namely electrodynamics; in Section 2.4, we make 
the connection to general relativity by (qualitatively) describing the major dif- 
ferences encountered in applying the described procedures to general relativity. 
In Section ||, we describe spin networks, a complete and sufficient basis for the 
states of loop quantization theory, as well as the actions of operators on this 
basis. 

Finally, in Sections ^ and we describe the results obtained in our research 
on the angle and volume operators, respectively. In particular, we find that 
the minimum observable angle for a given vertex is inversely proportional to 
the square root of its total spin — a fairly slow convergence to zero. We also 
calculate the "resolution" of angles around a vertex given its total spin, and 
show that the distribution of angles around a vertex with sufficient total spin 
approximates the classical distribution. The complexity of the volume operator 
hinders any serious analytical work; however, we do use numerical methods to 
show that the volume of a region seems to scale as the 3/2 power of the area of 
its bounding surface, which is also expected classically. 
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2 Theory of Loop Quantization 



2.1 Canonical Electrodynamics 

As stated in the introduction, we will begin by examining loop quantization 
theory in the context of classical electrodynamics. The fundamental equations 
of electrodynamics are, of course. Maxwell's equations: 

V-E = VxB-i|| = 

c at ,2.1) 

V • B = V X E+ = 

(note that here and throughout this section, we will be working with the source- 
free equations in Gaussian units.) The symmetry of these equations is striking; 
if we replace E with B and B with — E, the equations are unchanged. One 
might wonder, then, what the source of this symmetry might be. 

To answer this question, we need to return to one of the fundamental prop- 
erties of E and B: their definition as vectors. We can define any vector x in 
terms of its (Cartesian) coordinates (x^, x^, a;-^).!^ In this case, 

yL = x^ei+ x^e2 + x^ey, (2.2) 

To save effort, we will introduce the Einstein summation convention here: if an 
index is repeated in an expression, then we interpret this as summation over 
all possible values of that index. Using this convention, the previous equation 
becomes 

X = x'ei (2.3) 

Special relativity, however, tells us that we must also consider an event's 
position in time (not just in space.) To take this into account, we modify 
the conventional idea of a vector to create a four-vector; instead of having 
a component for each of the a;-, y-, and z-directions, a four-vector also has a 
component for the i-direction. We will denote these components 
x, = y, and = z. An event with spacetime coordinates {x^, x^,x^,x^) can 
then be written as 

X = x^e^ (2.4) 

In general, we will use Greek indices such as /z, i', . . . to denote components that 
must be summed over space and time, while Roman indices such as i,j, k, . . . 
will denote summation over space only. 

Special relativity tells us that if we switch between two frames travelling at 
a relative velocity P — v/c to each other in the x-direction, the components of 

^ Note that these superscripts are merely indices, not powers of x. 
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a spacetime event x will transform as 



(2.5) 



where the "barred" components are in the frame to which we are transforming. 
Using the summation convention outlined above, we can write (|]^) as a single 
equation: 



where A is the Lorentz transformation matrix 



A 



/ 7 

-7/3 


V 



-7/3 0^ 

7 

1 

1/ 



(2.6) 



(2.7) 



and AJ^ is the component of the /ith row and the vt\\ column of this matrix. 

If E and B are really vectors, then, they should transform according to 
equation (|]^). However, a Lorentz transformation on a given electromagnetic 
field gives the components of the transformed fields as 



Ex 



Ex 
Bx 



Ey 

B,, 



jiEy-pB,) 

-f{By + [3E,) 



E, 
B, 



j{E,+l3By) 

^{B,-pEy) 



(2.8) 



— not at all a simple Lorentz transformation like those mentioned above. 

The way that the components of E and B mix suggests that these two fields 
are really part of a larger entity. As it happens, they are components of a second- 
rank tensor. Such a tensor has sixteen components, which we can denote in an 
arrayiH 



^00 -foi -^02 -Fo3 

^ ^10 Ell Fi2 Fi3 

' F20 E21 F22 E23 

E30 E^i F32 F33 



(2.9) 



Theoretically, each one of these components could be different; however, such an 
object has 16 degrees of freedom, whereas we only have six to account for from E 
and B. As it happens, if we require F^^^ to be antisymmetric (i.e. F^^^ — — i^,y^), 

Note that, technically, this tensor should not be considered a m atrix, since it has two 
lowered indices; if we were to multiply it by a four- vector as in ( [2.6[ ), we would not obtain 
another vector j/". 
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such a tensor would only have six degrees of freedom: 






Fqi 


Fo2 


^03 







Fl2 


^13 


Fo2 


— Fi2 





F23 


Fo3 


~Fi3 








(2.10) 



How would the components of such a tensor transform under a Lorentz 
transformation? In analogy to the transformation of vectors in (2.6), we have 



Ffj.u - A'^A^Ff,^ 



(2.11) 



If we write out the transformations of each component of f^jy explicitly, we 
obtain 



Fpi — Fqi F02 — 7 (-^02 ~ PF12) F03 — 7 (^03 + PF^i) 

F23 — F23 F31 — 7 (i^3i — (3F03) F12 — 7 {F12 + /3F02) 



(2.12) 



The reader will note that these are exactly the transformations stated in (2.8), 
suggesting that we can write the electric and magnetic fields together in the 
field tensor F^^, as 



F„, 





-Ex 

-Ey 

-E,_ 



Ex 




E. 



V 



E, 



B, 
-Bx 



(2.13) 



We have thus succeeded in combining the electric and magnetic fields into a sin- 
gle tensor F^j_i, ; the peculiar transformations of E and B are merely consequences 
of the properties of this tensor. However, the symmetry in Maxwell's equations 
suggests that this is not the only way to put the components of E and B into a 
second-rank antisymmetric tensor; we could equally well have constructed the 
dual tensor G^": 





—Bx 
-B„ 



V 



Bx B, 
-E, 
E, 



-B, -E, 



V 



B 

E: 

-Ex 
Ex 



V 



(2.14) 



In terms of these new tensors. Maxwell's equations become radically simpler. 
They are now: 



d^F"'' = and d^G^"" = 0, 



(2.15) 



where we have introduced the operators 9^ — d/dx^ . That these equations are 
equivalent to Maxwell's equations is not entirely evident at first glance, so let 
us compute a couple of the equations (note that each of the equations above is 
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really four separate equations, one for each value of n) to show this equivalence. 
Consider the first equation when /i = 0; this becomes 



^ dE, ^ dEy ^ dE, 

dx dy dz 
= V-E = 0, (2.16) 

which is Gauss' Law. In the case where /i = 1, we have 

duF^" = daF^° + diF^^ + d2F^^ + dsF^^ 
IdE^ dB, dBy 
c dt dy dz 

=0 (2.17) 

c dt 

which is the x-component of Ampere's law. The /x = 2 and /i = 3 cases of this 
equation give the y and z components of Ampere's law, and the tensor equation 
for G'"' gives the two remaining Maxwell's equations 



2.2 The Vector Potential, Holonomies &6 the Aharonov- 
Bohm effect 

To quantize a given theory, we need to decide on a particular set of configuration 
variables which describe the system. In the simplest example, the quantum the- 
ory of a free particle uses the particle's position x and momentum p to describe 
the motion of the particle. In general, we want to find some set of generalized 
"positions" and "momenta" that can describe all the states of the system; how- 
ever, doing so is not always easy, since a "natural" set of configuration variables 
may be related by some set of complicated constraints. A general method for 
finding the constraints on an arbitrary set of configuration variables was pio- 
neered by Bergmann and Dirac, and is presented in Appendix along with its 
application to electrodynamics. Here, however, we only need to note the final 
results: the canonical configuration variable for electrodynamics is the vector 
potential A, given by B = V x A. We will also make use of the four-vector 
potential, given by 

A^ = {V,A,,Ay,A,) (2.18) 

where V is the electric potential. The field tensor i^^^ can be written quite 
elegantly in terms of this four-vector potential: 

F^, = d^,A, - d,A^. (2.19) 

* This method can be elegantly extended to include field sources (i.e. charges and currents), 
but this case will not concern us here. The interested reader is referred to Grifliiths [h| (from 
which most of this derivation is taken) for further information. 



7 



To successfully quantize this theory, we must also have a conjugate momentum 
for the chosen configuration variable. Using Bergmann-Dirac canonical analysis, 
we can show that the "momentum" of the vector potential Aa is the electric field 
E''!^ Our two configuration variables are the vector potential and the electric 
field. 

One might ask the question, however, whether the vector potential is nec- 
essarily the best choice for a configuration variable, especially since (in the 
classical theory) only B = V x A is measurable. Furthermore, Aa has the dis- 
tinct disadvantage of possessing gauge freedom: given a vector potential Aa , we 
can perform a transformation of the form 

. . d\ 

without affecting the physical state. This freedom is worse than it appears: 
suppose we have a physical system which we consider under two gauges Ai and 
A2, the first of which is a function of time and the second of which is not. If we 
let the system evolve, the two gauges will generally predict different values of 
Aa after a given amount of time; hence, our theory does not uniquely predict 
the time-evolution of the configuration of the system. If we insist on treating 
Aa as a fundamental variable, then, we must abandon determinism — a step 
that we are loath to take. Further, this interpretation does not prohibit causal 
influences from travelling faster than light; in theory, a change in the field can 
cause an instantaneous change arbitrarily far away. Thus, interpreting Aa as a 
fundamental field forces us to eliminate determinism and causality, two of the 
most basic tenets of physics. 

In the face of this, most people would interpret B = V x A (and E) as 
physical fields; this is the standard classical interpretation. This interpretation 
has the advantage of being deterministic (since gauge transformations do not 
change B) and limiting causal influences to the speed of light. Indeed, if we 
are defining our fields on a "well-behaved" chunk of space, or only examining 
classical phenomena, then the choice of B as the fundamental variable is natural, 
and explains everything that needs explaining.^ 

As the reader might have guessed from the above caveat, however, there are 
quantum effects that seem to contradict the primacy of B. The best- known 
of these is the Aharonov-Bohm effect. In this experiment, a beam of charged 
particles with charge e is sent through a double-slit interference experiment 
(see Figure ^) Between the slits we place an impenetrable solenoid (i.e. the 
particles' wave functions are zero inside the solenoid.) 

Theory predicts that a particle travelling along a given path 7 will change 

^ Note that in this equation, the components of E have a superscript instead of a subscript; 
this is due to the difference between covariant and contravariant tensors. In this case, however, 
raising or lowering an index simply introduces a minus sign if and only if the index is (i.e. 
x'^ = —xo, but x^ = XI.) 

^ The reader might wonder what we might possibly mean by "well-behaved"; this will 
become evident shortly. 
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B = 




Screen 

Figure 1: The Aharonov-Bohm experiment. Particles are sent through a 
double-sht apparatus and create an interference pattern on the screen. The 
presence of the solenoid creates a phase difference between the two paths 71 
and 72, which affects the interference pattern. Note that although the magnetic 
field is zero in the region where the particle travels, the solenoid still has an 
effect on the particles. 
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its phase factor: 



\(f>{t)) = exp 



-iEt 



+ 



le 
he 



(2.20) 



where E is the energy of the particles. The first part of the phase factor corre- 
sponds to the normal time-evolution of the state; the second part arises solely 
from the vector potential due to the solenoid (even though the particle is mov- 
ing in a region of no magnetic field.) This phase factor will have observable 
consequences in the interference pattern on the screen behind the slits. For 
example, the phase difference between the two paths at the centre of the screen 
(i.e. equal time for both paths) will be given by 



= ^ / A,dx'-^ 



71 



72 



e 

he 



Ai dx' 



(2.21) 



where 7 is the closed loop created by travelling along 71 , then along 72 in reverse. 
The integrand of this equation can be expressed more familiarly as A • dx; we 
can thus apply Stokes' Law to this to give 



A(/) = — / V X A • dS 



e 

he 
e$ 
he 



(2.22) 



where 7 is the border of the new surface S, and $ = / B • dS is the magnetic 
flux through this surface (which is equal to the flux through the solenoid, since 
B = outside it.) 

This experiment has been performed (the experimenters measured shifts in 
the interference pattern as $ was increased) ; the results confirm the theoretical 
prediction, that the interference pattern does indeed depend on J A • dx. This 
result seems to catch us in a contradiction. On the one hand, choosing the 
vector potential A^ as a fundamental configuration variable forces us to aban- 
don determinism and causality; on the other hand, the Aharonov-Bohm effect 
suggests that the vector potential holds some physical importance. Without a 
third option, we are forced to choose the former — while a non-deterministic 
theory is not very nice, one that doesn't correspond to the physical world is 
worthless. 

Fortunately, a third option does exist: holonomies. Recall that the observ- 
able result of the Aharonov-Bohm effect is not depend ent o n the vector potential 
per se, but instead on its path integral as defined in (2.21). Thus, given a path 
in space 7, we define the holonomy of 7 to be 



exp 



exp 



W / A,ijit))Y{t)dt 



ig I Ai dx^ 

7 



(2.23) 
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where 17 is a coupling constant to make the argument of the exponential dimen- 
sionless. The reader will note that if we set g = e/hc, this factor reduces to that 
given in the second term of (2.20). 

We can thus escape our earlier dilemma by regarding the holonomies as 
the fundamental configuration variables. The consequences of doing this are 
discussed at length in an excellent article by Belot He notes that the 

holonomies, like B, are gauge-invariant; however, they possess another impor- 
tant, related quality which the magnetic field does not. In a chunk of "well- 
behaved" space, two fields correspond to the same B field if and only if they 
differ by one of the usual gauge transformations (i.e. a term of the form dX/dx"" 
for some A.) This is because normal Euclidean space is simply connected: any 
loop in space (i.e. a path through space that begins and ends at the same point) 
can be contracted continuously to a point. However, if a space is not simply 
connected, then there are non-gauge transformations that are interpreted as 
identical B fields. The holonomies, on the other hand, do not have this prob- 
lem; this is thus a point in their favour.^ 

However, using the holonomies as the fundamental variables of the theory 
violates locality in an even more fundamental way than our first interpretation. 
Instead of considering the field Aa at a collection of points, we must now consider 
the values of the field over a region with actual spatial extent; this means that 
if we don't know the field everywhere^ we cannot determine all of the possible 
holonomies. The philosophical issues this fact raises are fascinating, but are 
unfortunately not the primary concern of this paper; the interested reader is 
encouraged to consult Belot's paper for more information on this subject. 



2.3 Holonomies and Reconstruction Theorems 

We have shown that the construction of a quantum theory from the canonical 
variables Aa and i?" depends in some way on holonomies; we will now make 
this notion somewhat more formal. We define a path 7 on a space E to be a 
continuous map from the closed unit interval [0, 1] to S. In the case of Euclidean 
space we can think of E as equivalent to R"^; however, E may equally well be a 
"curved" manifold, as in general relativity. If the property 7(0) = 7(1) holds, 
then we say that 7 is a closed loop. 

Given a (not necessarily closed) path 7 and a vector potential Aa^ we can 
then define a holonomy U.y of the path 7: 



— exp 



exp 



ig I v4,(7(t))7°(i)d< 




ig I Aa dx"' 
7 



(2.24) 



where g is merely a coupling constant introduced to make the argument of the 

^The reader who might be wondering if non-simply-connected spaces actually correspond 
to any physical situations is reminded that the Aharonov-Bohm effect takes place in such a 
space — a loop around the solenoid cannot be contracted "through" the solenoid. 
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exponential dimensionless. In other words, the holonomy of a given path is 
merely (in this simple case) the exponential of i times the path integral of the 
vector potential along the path in question. 

An important feature of closed-loop holonomies is that it is generally possible 
to reconstruct the vector potential from them; in other words, the representation 
of the vector potential in terms of its holonomies and in terms of the field itself 
are completely equivalent. The conditions under which this holds are the subject 
of mathematical reconstruction theorems. We will not go into detail on how the 
existence of a vector potential corresponding to a given set of loops is proven, 
nor the actual process of the reconstruction itself; for more details, the interested 
reader is referred to Barrett |3| or Giles jij. 



2.4 Canonical variables &: holonomies in general relativity 



Before we examine the nature of the spin networks mentioned in Section 2.3 
we must discuss the difficulties in applying this technique to Einstein's the- 
ory of general relativity. As originally conceived by Einstein, general relativity 
consisted of a spacetime manifold S, at each point of which a metric g^i, is 
defined. Early attempts to quantize general relativity used the metric as one of 
the parameters of the configuration space. Unfortunately, this approach failed, 
mainly because the constraints generated by this approach (as in Appendix ^ 
were intractable. 

There the matter rested until the 1986, when Ashtekar ^ and Sen [Q 
reformulated general relativity in terms of the so-called "new variables" . In this 
formulation, the metric was no longer the fundamental variable of the theory. 
Instead, the two variables used in this new formulation are the connection A]^{x) 
(akin to the vector potential in electrodynamics) and the inverse densitized triad 
Ef{x) (akin to the electric field). As our notation would suggest, the connection 
is viewed as the configuration variable (akin to the vector potential), while the 
inverse densitized triad is its canonically conjugate momentum. These quantities 
have "clean" Poisson brackets, as we require: 

{Et{x),E]{y)} = (2.25) 

{Aj.(x),^^,(y)} = (2.26) 

{A^Sx),E]{y)] ^ (3GSiSl6'{x-y) (2.27) 

For those readers familiar with general relativity, we can also express these 
variables in terms of more familiar quantities. Ef is related to the metric of 
constant-time surfaces Qat by 

Iqlq^" = E^El (2.28) 

where qab is the spatial part of the metric and \q\ is the determinant of qat- We 
can also express these in terms of the triad from the ADM formalism: 

= VWK (2.29) 
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The connection Al^{x) can be expressed as 

Al^Tl+f3Kl (2.30) 

where FJ^ is the connection associated with the (undensitized) triad (defined by 
9[a^b] ~ ^la^b]j)i cxtrinsic curvature of the 3-surface, and P is an as- 

yet undetermined parameter known as the Immirzi parameter, first put forward 
in §. 

While there is a resemblance between electrodynamics and general relativity 
in this regard, there are also differences (of course) between these two theories. 
The reader will no doubt have noticed by now that our relativistic variables 
bear an extra index (i.e. . . .). This index corresponds to the internal space 
of the configuration variables. Each point on E has an internal space associated 
with it; a vector based at a given point can be thought of as existing in this 
space. Since two vectors at two different points exist in different spaces, they 
cannot immediately be compared. We need to define some self-consistent way 
of comparing two objects at two different points; this, as it turns out, is done 
by the connection A^. 

Suppose, then, we take a vector at some point xq on the manifold, and at- 
tempt to move it around the loop (using the connection to translate from 
the neighbourhood of each point to the manifold.) This movement will be per- 
formed in such a way that as the vector passes through a small neighbourhood 
of each point, it will point in the same direction. This process is called par- 
allel transport, since the vector stays parallel to itself along each infinitesimal 
segment. 

Transporting a vector around a loop in this way does not always yield a 
vector that is parallel to the original vector — even if we keep the vector parallel 
to itself along each "step of the way." In the special case of flat space, the 
parallel-transported vector is equivalent to the original vector. However, if space 
has some amount of curvature, we will not necessarily get the same vector back 
after we transport it around a loop. (For an example of this, see Figure ||, in 
which a vector is parallel-transported around two different loops on a sphere.) 
With this in mind, we associate a three-dimensional rotation with each path 
in E; these rotations will be the holonomies of general relativity. In three 
dimensions, we can represent these rotations by members of the group 5*0(3), 
the group of all orthogonal 3x3 matrices with determinant 1. However, this 
group has the same structure as (i.e. is isomorphic to) SU{2), the group of all 
2x2 unitary matrices with determinant 1; since these matrices are simpler, 
general relativity is more commonly though of in terms of SU{2) than in terms 
of 50(3). 

This representation of the holonomies in terms of matrices has two impor- 
tant consequences that complicate the simple derivation we used for canonical 
electrodynamics. The first complication arises from the structures of the groups 
used in each theory. In the case of canonical electrodynamics, the holonomies 
were complex numbers of unit modulus; these can be thought of, in the interest 
of drawing parallels, as members of U{1) (the set of all unitary 1x1 matrices.) 



13 




(a) (b) 

Figure 2: Parallel transport of vectors on a non-flat surface. In (a), the vector 
is transported from the north pole to the equator, 90 degrees along the equator, 
and back to the north pole; the resulting vector is rotated by 90 degrees. In (b), 
the vector is transported from the north pole to the south pole along a line of 
longitude, and back to the north pole along another line of longitude 90 degrees 
from the first; the resulting vector is rotated by 180 degrees. Note that the 
vector has been transported "parallel to itself" at each point; it is the curvature 
of the surface itself that causes these rotations. 



U{1) happens to be a simple group, in that any two elements a,b d U{1) com- 
mute (i.e. ab = ba.) However, SU{2) is not so simple; in general, elements of 
SU{2) do not commute (i.e. ab ^ ba.)'^ Beca use of this, it turns out that we 
can no longer use the simple equation in ( 2.23 ) to find the holonomy. Instead of 
simply taking the exponential of a path integral, we must use the path-ordered 
exponential defined as 



Uj[A] = -Pexp 



Aah{t))rit)dt 



(2.31) 



1 + 



dti 



1 

dt2' 



dtr. 



xAaSl{tn)) ■ ■ ■ Aadl{tl))T"{tn) ■ ■ ■ i"' (tl) (2-32) 



Readers familiar with time-dependent perturbation theory may recognize this 
expression as being similar to the time-ordered exponential (also known as the 
Dyson series.) 

The second important complication introduced by the matrix representation 
arises when we attempt to account for gauge transformations. Under gauge 
transformations, the holonomies transform as 



7 



^=5~'(7(l))-f/7- 5(7(0)) 



(2.33) 



* The reader with a modicum of background in group theory will recognize that U{1) is 
Abelian, while SU{2) is non-Abelian, 
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where 5 is a smooth SU{2)-vahicd function on E. This kind of thing just won't 
do; we want a gauge-invariant quantity. To turn this quantity into such a quan- 
tity, we restrict ourselves to closed loops, and take the trace of the holonomy, 
also known as the Wilson loop: 

= Ti{U-y) (2.34) 

This quantity is invariant under gauge transformations: 

= Tr(.g-i(7(l))-C/v .9(7(0))) 
= Tr(ff(7(0)) .,9-^7(1)) -C/^) 

= Tt{U^) = T-,, (2.35) 

where we have used the cyclic properties of the trace and the fact that 7(0) = 
7(1)- 
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3 Spin Networks 



The construction of a quantum theory from closed-loop holonomies is known 
as the loop quantization technique. However, there is a problem with loop 
quantization: the loop basis is overcomplete. In other words, certain linear 
combinations of closed-loop holonomies are always equal to each other. This 
means that if we construct a theory out of these loops, we will end up with 
linearly dependent basis states (since we will not have taken these equalities 
into account.) 

In attempting to resolve this difficulty, Rovelli and Sniolin |9j discovered 
that it could be circumvented entirely by the use of spin networks. To construct 
such a network, we must first make another definition: we define a graph T as 
a finite collection {71,72, . . . ,7™} of smooth curves in E that intersect only at 
their ends (if at all.) With each curve, we associate its holonomy, i.e. 



U^^ (A) = V exp 



A„ Ax'' 



7i 



(3.1) 



Our quantum states will then be functions of these holonomies: 

*(A) -/([/i,C/2,... ,C/„) (3.2) 

Since these networks play a central role in our theory, we will take some time 
to examine them now. Much of the following exposition is taken from Major 

m. 



3.1 Definition & Properties 

Spin networks were first invented by Penrose, ||ll| and essentially consist of 
graphs with weighted edges. Their fundamental units are strands, which are 
defined to be equivalent to 2 x 2 matrices: 

B 



A~ / e = ^e^"" = \j e = i^AB = r\ (3.3) 
A ^ ^ 

Indeed, any 2x2 matrix can be "placed" into a spin network, in the form of a 
"tagged" fine: 



B 



^ (3.4) 

A 

To join two strands together, we multiply two matrices with a common index: 

/ ' 

6^ ■ ■ ecD ■ 5E = B^cJ = ^ = 5| (3.5) 
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Using these rules, we can show that these strands behave as would thin 
strings in a plane, with one exception: the equation S^Sj^ecD — —^ab leads us 
to the odd strand equation 

%/ = ^A^coS'i = Iba = -~€AB = - r\ (3-6) 

r\ A B 

A B 

To resolve this, we associate a negative sign with each strand crossing. With 
this modification, the strands of a spin network behave exactly as would thin 
(i.e. untwistable) threads in a plane. 

The mathematical definition of the strands in terms of matrices gives rise to 
an interesting identity, known as the "binor identity" or the "skein relation." 
We recall the matrix identity 

e^ce^^ = 5^5^ - 6^6^ (3.7) 

When translated into strands, this identity becomes 

B D B D 

= (3.8) 



The fact that a given set of strands may satisfy one or more linear relations 
(as in the binor identity above) leads us to construct a basis in which such 
linear relations do not exist. This is accomplished by the use of antisymmetric 
combinations of multiple strands. The simplest example of such a combination 
is the 2-strand: 




I 1 

f=2 



(3.9) 



In general, we construct an n-strand by summing over all permutations of the 
lines, changing the sign for each strand crossing: 



1 = ^ E(-i)''^' tzn (3.10) 

I ■ o-es„ I I 

where a stands for a permutation of the lines.|^ 

^ The exact correspondence between elements of Sn and "permutations of the lines" can 
be established as follows: represent an element of Sn by its permutation on n letters. Write 
the letters 12 ... n in a row, and write the permuted letters a(l)cr{2) . . . (j{n) just below them. 
The lines connecting the same elements in the top and bottom rows correspond to the weaving 
of the strands. 
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These networks would be boring if they only consisted of simple n-strands. 
However, we can define a trivalent vertex of these edges by splitting three n- 
strands and reassembling them as follows: 

\ / 
\. ./ 

¥ 

c 

I 

where the labels i, j, and k are "internal" labels, and are related to the "exter- 
nal" labels a, b, and c by the relations 

i = i(a + c-6) 

j = 1(6 + c- a) 

k = ^{a + b — c) 

The dotted lines in the above diagrams indicate the boundaries between the 
internal structure of the vertices and the external structure (i.e. the strands 
themselves.) For these internal edges (and hence the vertex itself) to exist, a, 
b, and c must satisfy two conditions. First, the sum a + b + c must be even; 
second, a, b, and c must satisfy the triangle inequalities: 

a + b>c b + c>a c + a>b 




3.2 Intertwiner Bases &; Recoupling Theory 

Even with the existence of trivalent vertices, life is rather boring; we ask our- 
selves, then, whether it is possible to construct vertices with arbitrary valence. 
We could try using the same technique we used in constructing the trivalent 
vertices, by "splitting" the incoming strands and reassembling them in all com- 
binations: 




However, we quickly notice that the internal edges are underdctermined: there 
are six internal edges i,j,k,l,m, and n, and only four external edges to de- 
termine them. Indeed, for a general n-valent vertex, there will be n{n — l)/2 
internal edges to be determined by only n parameters. 
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While we could still proceed with our analysis using all of the extra internal 
edges, it will be to our advantage to introduce a different way of denoting 
the internal edges of a vertex. Since internal edges of a trivalent vertex are 
completely determined, we will attempt to build a four-valent vertex out of 
trivalent vertices instead: 

^1. 



Instead of having to deal with a mishmash of internal edges, we will now only 
only have to deal with the one internal edge i. Indeed, from here on the phrase 
"internal edges" will be used exclusively to refer to edges such as i, and not to 
the "split strands" used to construct the trivalent vertices. 

However, this construction is useless if we cannot actually find an edge i that 
allows the internal trivalent vertices to be constructed. The first requirement to 
deal with is the "evenness" requirement. For the trivalent vertices to exist, the 
restrictions on the triples (o, b, i) and (c, d, i) are: 

a + b = i (mod 2) and c + d = i (mod 2), (3.11) 

which together imply that the sum a + h + c + d must be even. Similarly, 
from the triangle inequalities a + b > i and i + c > d together imply that 
a + b + c > d; we can similarly construct three other inequalities, creating a 
total of four "generalized" triangle inequalities: 

a + b + c> d 

a + b + d> c , ^ 

3.12 

a+c+d>b 
b + c + d> a 

Alternately, these four inequalities can be summarized in one statement: 

a + b + c + d> 2x for X e {a,b,c,d} (3.13) 

There is still an ambiguity in the internal edge, however; a bit of thought 
reveals that the value of i can lie in the range 

max{|a - b\, \c- d\} <i< min{a + b,c + d} (3-14) 

For example, in the case where a = b = c = d, i can take on any (even) value 
between and 2a. We can think of these sets of vertices as state vectors, which 
span the space of all possible four-valent vertices with "external" edges a, b, c, 
and d. 

We can define the inner product between two vertex states by joining to- 
gether corresponding external edges. For example, if we take the inner product 
between a state labelled by i and a state labelled by i', we find that it is given 
by the value of the following network: 
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U U L U 



If we apply the "bubble" identity ( |B.5| ) to right-hand bubble in this diagram, 
we find that the value of this diagram is unless i — i': 

a b c d ~ Ow 7 a b i 



9{i, a, b)0{i, c, o?) 



(3.15) 



Thus, these basis states are orthogonal.0 

The reader may have wondered, back when we were defining the four-valent 
vertex in terms of internal trivalent vertices, why we chose to do so in this 
manner. After all, we could have equally well joined together the vertices a and 
d in one internal vertex and the vertices b and c in the other: 



If we had defined the internal edge in this way, we would have obtained a differ- 
ent set of state vectors; however, it is not immediately evident that these vectors 
will span the same space as our original set. That these new vectors actually 
do span the same space is a result of the recoupling theorem (see Appendix 



d 

/ \ 



Thus, given a vertex with an internal decomposition in either form, we can 
express it as a linear combination of vertices with an internal decomposition of 
the other form. These internal decompositions are called intertwiners; what we 
have just shown is that the set of all intertwiners of either form is a basis for the 
space of all intertwiners, and that the recoupling theorem is a transformation 
between these bases. Moreover, this transformation happens to have a nice 
property. In jl2j, De Pietri and Rovelli show that the recoupling theorem is 
a unitary transformation, i.e. vector norms are preserved under it. We can 
thus switch back and forth between bases fairly simply; this ability will become 



"They are certainly not normalized; however, this fact will have no effect on the properties 
we examine. 
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Figure 3: The "branch transformation", in which an internal branch is added 
to the intertwiner. Internal edges are marked with is and js, while external 
edges are marked with as. The gray outlines indicate the internal edge to which 
the recoupling theorem is applied in the next step. 



invaluable to us when we examine the actions of the angle and volume operators 



in Sections 3.5.3 and 3.5.4, respectively. 



We have now defined what we mean by an 4-valent vertex; however, our 
goal was to create a vertex of arbitrary valence. It is fairly obvious how this 
can be done: for each new edge we wish to add, we add a trivalent vertex to 
the intertwiner. Thus, a 5-valent vertex contains three trivalent vertices in its 
intertwiner, a 6-valent vertex contains four trivalent vertices, and (in general) 
an n-valent vertex contains n — 2 trivalent vertices in its intertwiner. It can also 
be shown that given a set of incoming edges with labels {ai, 02, . . . , a„}, it is 
possible to create an n-valent vertex out of them if and only if the following 
conditions are satisfied: 



E 



Qi IS even 



(3.17) 



^ai> 2aj for 1 < j < n (3.18) 

i=l 

The first of these is the usual "evenness" requirement, while the second is a 



generalized triangle inequality similar to ( 3.13 ). 

As might be expected, these higher-valence vertices also have intertwiner 
bases (although these bases are of much higher dimension.) However, we can 
again use the recoupling theorem to shift between these bases. There are two 
general types of these transformations which will be very useful to us later on. 
The first of these is illustrated in Figure ||, and allows us to transform between 
bases of different "shapes." By creating (or, in reverse, removing) branches 
from an arbitrary internal decomposition, we can transform an intertwiner of 
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Figure 4: Two useful shapes for intertwiner bases. The first, in (a), is the 
"snowflake" basis. Its essential feature is that the edges of a given n-valent 
vertex are partitioned into three groups (denoted {xj}, {yi}, and {zj} in the 
diagram); each of the groups is then "collected" into a single edge (labelled ni, 
712, and Ur); and finally, these three collecting edges form a trivalent vertex at 
the core of the vertex. The second basis, in (b), is the "comb" basis; in it, all 
the edges are lined up along a central "spine" of intertwiner edges (labelled 12, 
13, ii . . . in the diagram.) 



arbitrary shape into any other shape we wish. Two such useful shapes are the 
"snowflake" and "comb" bases, depicted in Figure ^. 

The second such transformation is illustrated in Figure ^. This transforma- 
tion allows us to switch the order of two adjacent edges along a branch of the 
intertwiner, while leaving the other branches in their original order. Thus, any 
order of edges is possible along a given intertwiner shape; this fact will become 



vitally important to us in our discussion of the volume operator in Section 3.5.4 



3.3 Spin Networks & Angular Momentum Coupling 



The trivalent vertex rules introduced in Section 3.1, and their generalized ver- 
sions m (p7|) and (|3l8| ), are not entirely transparent. However, a simple 
analogy exists which may make these rules more intuitive to the reader; there 
is a striking parallel between the coupling of angular momenta in quantum me- 
chanics and spin networks p| 

To draw this parallel, let us relabel all the edges of a given trivalent vertex 
with half their original values. In other words, an edge with colour 1 becomes an 
edge with colour ^; an edge with colour 2 becomes an edge with colour 1; and 



in general, an edge with colour i is reassigned a new colour i' - 

^^Indeed, this similarity is the very source of the name "spin" networks. 



i/2. In terms 
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J2 fc2 



fli 02 fc2 



Figure 5: The "swap transformation", in which two edges along an internal 
spine are s witched. Note that in the third step we have applied the lambda 
move (B^) to the vertex indicated in the previous step. 



of these half-integer edge colours, the constraints that the edges must satisfy 
become: 



a' + b' > c' b' + c' > a' c' + a' > b' (as before) 



(3.19) 



b' + c' is an integer 



(3.20) 



A it happens, these rules are exactly the same as those for the coupling of 
particles with angular momentum. If we wish to couple together a spin-a' and 
a spin-&' particle to form a spin-c' particle, the three spins a', 6', and c' must 



satisfy the triangle inequalities, as in (3.19); as well, the sum of the three spins 
must be an integer, as in ( 3.20| ). 

As a concrete example of this parallel, consider a trivalent vertex with an 
edge with colour 1, an edge with colour 2, and the third edge colour m as yet 
undetermined. According to the vertex rules, we must have m = 1 or m = 3; any 
value less than 1 or greater than 3 would violate the triangle inequalities, and 
a value of 2 would violate the evenness requirement. The angular momentum 
coupling scheme which correspond to this vertex would be the coupling of a 
spin-i particle to a spin-1 particle; according to the rules of such couplings, the 
resulting spin must be either ^ or ^. 

We can draw a similar parallel between 4-valent vertices and three-particle 
coupling schemes. Consider again the following vertex: 



Suppose a = 2, & = 3, and c = 3. According to the vertex rules, this means 
that d must have a value of 0, 2, 4, 6, or 8. Similarly, if we couple together 
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one spin-1 particle and two spin-| particles, the resulting total spin can be 0, 
1, 2, 3, or 4. In fact, this rule holds for a general vertex with valence n, with 
incoming spins ai, 02, • • • , a„; such a vertex can be translated into a coupling 
scheme of rt — 1 spins a'i,a'2, . . . , a'n-i to yield a total spin of , where a[ — ai /2 
for j = 1, 2, . . . ,n. 

The reader may be asking at this point whether any analogy for the internal 
edges of a higher- valence vertex exists in angular momentum coupling schemes. 
In fact, such an analogy does exist. In the above diagram, we can view the total 
coupling as happening in two stages: first particle a and particle b couple to 
form a particle i, then particle i and particle c couple to form particle d. Of 
course, the coupling does not necessarily have to take place in this order. For 
instance, we could equally well have coupled particles b and c together first, to 
form a particle j, and then couple particle j to particle a to form particle d; 
this is an equally valid way to couple the three particles together. Similarly, we 
can represent a vertex with edges a, 6, c, and d in a second form: 




The parallel between the intertwiner of this vertex and this second angular 
momentum coupling scheme is evident. 

Finally, we note the following: In angular momentum coupling, we translate 
between these two coupling schemes (i.e. particle i vs. particle j) using the 
Wigner 6-j symbols; in spin networks, we use the recoupling theorem and the 
Kauffman-Lins 6-j symbols to translate between intertwiners. Given the above 
parallels, it should not surprise us to learn that the Wigner and the Kauffman- 
Lins 6-J symbols are related; for more details, see Appendix 

3.4 Spin Networks & Holonomies 

The reader may be wondering at this point what, exactly, the connection be- 
tween the spin networks we have just described and the holonomies described 
in Section || is. To answer this question, consider a simple network in space, 
shown in Figure ^, along with three paths (or sets of paths) in this network. To 
label these paths, we must introduce a few concepts. We define the composition 
of two paths to be the path obtained by traversing the two paths in order; for 
example, the path a o 7 is obtained by going around a, then going along 7. We 
also introduce the concept of inverses of a path, which we define (intuitively) 
to be the same path traversed in the opposite direction; for example, the path 
a o 7 o 7~^ is obtained by traversing a forwards, then 7 forwards, then 7 back- 
wards (i.e. retracing our steps along 7). Using this notation, the three paths 
can be labelled aU (3, ao j o (3 o and a o 7 o o 7"^. 

One might think that the state vectors corresponding to these three loops 
would be independent; unfortunately, this is not the case. Recall that each path 
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a 






a o ^ o p o J 1 




Figure 6: A simple network in space, and three paths (or sets of paths) within 
it. 



in the manifold E has a matrix from the group SU{2) associated with it. Given 
any two matrices A,Be SU{2), their traces satisfy the following identity: 



Tr(A)Tr(B) - Tt{AB) - Tt{AB-^) 







(3.21) 



Since this linear combination is equal to zero, the quantum states corresponding 
to each one of these terms satisfy an analogous identity: 



(a U /3| — (a o 7 o o 7 ^\ — (^aojop o 7 -'-|=0. 



(3.22) 



Hence, these three states are linearly dependent. 

For a fairly simple network (such as the one just described), it might be 
possible to keep track of all the identities satisfied by loops in the network. 
However, for a more complicated network the number of these identities sky- 
rockets; it would be nice, then, if we could find a way to avoid these identities 
altogether. 

As the reader may have guessed by now, spin networks allow us to do this. 
Consider the states 



(q;U/3| = (a 070/307 "'"1 + ^01070/3 -'-07 
(^a o J o p o 7^""" I — o 7 o p~^ o I . 



and 



(3.23) 



These two states are independent of each other, and (just as importantly) span 
the entire space of loops on this network — any loop or set of loops on this 
network can be written as a linear combination of these states (using the identity 
in 3.22 if necessary.) However, the second one of these states can be rewritten 
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in a more familiar form: 




Hence, given a network in E, the introduction of spins on the edges ehminates the 
identities that might otherwise be present. In this way, spin networks represent 
a basis for the states of loop quantum gravity. 
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3.5 Operators on Spin Networks 

We have now defined our quantum states in terms of spin networks. A successful 
theory of quantum gravity should be able to make predictions about "spatial 
observables" such as area, volume, and angle — geometric quantities that we 
would normally be able to measure on a classical manifold. As in any quantum 
theory, this is done by Hermitian operators which act on the spin networks and 
return eigenvalues corresponding to observable quantities. However, it is not at 
all clear what these operators should be. 

We will describe the actions of four operators — angular momentum, area, 
angle, and volume — in varying degrees of detail. In general, however, the 
construction of quantum operators from first principles is a complex process; 
hence, we will not try to do so in this work. The reader wanting further details 
as to the construction of these operators is referred to the derivations in (for 
angular momentum), jl3| Q (for area), [|l5j (for angle), and [|l2j (for volume). 



3.5.1 Angular Momentum Operators 

One of our primary concerns in this thesis will be to describe the angles as- 
sociated with a given vertex; to define angles, however, we need to be able to 
associate a direction (i.e. a vector) with an edge or a group of edges. In Pen- 
rose's original formulation of spin networks, the labels on the edges actually 



refer to angular momenta "carried" by the edges (see Section 3.3); hence, one of 
the most natural directions to associate with an edge is its angular momentum 
vector. 

These angular momenta can be measured by the angular momentum opera- 
tors, which act upon the edges. The angular momentum operators are expressed 
in terms of the Pauli spin matrices (7^: 

-(;;)^ '-("o')' -(J-i) <-^' 

The operators themselves are then defined in terms of these matrices "grasping" 
the edges: 



-2- (3.25) 



where the small square at the end indicates that this operator will "grasp" 
the edge. If the strand in question happens to be a multiple strand, then this 
grasping will apply to each strand in it; hence, we can show that a 2-edge 
grasping an m-edge is equivalent to m times a trivalent vertex with edges m, 
m, and 2: 



-2-a in = m 



m 

2-| (3.26) 
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Using this, we can define angular momentum operators Ji = |(fi which act 
on the edges; these operators have the usual properties of angular momentum 
operators. As we might expect, we can also construct the operator. This is 
equal to the sum of the squares of the three Ji operators, i.e. = Jl + J^+J"^. 
We can use diagrammatics and graspings to find the eigenvalues of this operator; 
however, it is much simpler to note that the Pauli matrices have the following 
identity, which is easily expressible in terms of spin networks: 




-X]- 



A C 

(3.27) 



Hence, instead of expressing in terms of the sum of three operators, we can 
express it very simply in diagrammatic form: 



-4 



(3.28) 

where we now grasp with two "hands" on the edge in question. Using the 



"bubble" identity (B.5), this becomes 



J2 



-m 2 m 



, (-l)™6l(77i,m,2) 



1 



Note that the factor of in the firs t step comes from the two graspings of the 
2-edge. If we apply the identity (B.3) to this equation, we get 



J m — —m 

I 2 m + 1 



(-1)™ (-l)"+i(m + 2)(m + 1) 



2m 



= -h 



2m(m + 2) 



(3.29) 



If we make the identification j = m/2, then the eigenvalues simplify to the easily 
recognizable form h^j{j + 1).^ This operator (and its eigenvalues) will become 
important in our discussion of the area and angle operators. 
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Figure 7: A surface E whose area we wish to find. The contribution to the 
area from the portion of the surface shown would be the sum of the contribution 
from each of the two edges ii and 12 intersecting this portion. 



3.5.2 The Area Operator 

In the spin-network version of quantum gravity, the area operator is the simplest 
to define. If we have a closed surface S which intersects our spin network in 
N points, its area is can be decomposed into a contribution from each of these 
intersections: 

N 



(3.30) 



i=l 

As it turns out, it is much easier to derive the action of the area-squared 
operator Af than to directly find the action of the area operator; this action 
turns out to be the double grasping of a two-edge times the square of the Planck 
length Iq, defined in (Pl|): 

if = ii (3.31) 

The reader will notice that we have already done the work to find the eigenvalues 



of this operator in the previous section; the result, from (3.29), is then 

Af i = ll L 1 I (3.32) 
The total area of E is then given by 



^q^^oEV^^^^^T^ (3.33) 

i 

where the summation of i is over all edges that intersect E, and mj is the edge 
label of the edge i. 

The actual eigenvalues of the operators Ji will not concern us here; however, it is worth 
noting that doing so requires (loosely speaking) giving orientations to the edges. Since this is 
not required in the spin-network formulation of quantum gravity, we will not go into it here. 

The reader will notice that the sign of this expression is incorrect. This stems from a 
difference between the properties of the "angular momentum" spin networks originally studied 
by Penrose and those used for quantum gravity. 
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Figure 8: Partition of the edges surrounding a vertex into three groups 5*1, S2, 
and 5*3. (The vertex itself is inside the sphere.) The angle operator will return 
the angle between the patches Si and 82- 



3.5.3 The Angle Operator 

The angle operator acts on a vertex as follows: we can imagine a surface encom- 
passing a vertex of arbitrary valence. We select two "patches" on the surface 
Si and S2, and call the remaining area of the surface Sr (as in Figure This 
defines a partition of the edges into three groups corresponding to the three 
patches Si, S2, and Sr- The intcrtwiner basis for this vertex is decomposed in 
such a way that all edges belonging to one group are "collected" into one internal 
edge; the three internal edges created in this way then join at a trivalent "core;" 
see Figure For future use, we will also define the quantities si, S2, and Sr 
as the sum of the labels of all the edges traversing Si, S2, and Sr, respectively; 
similarly, we will define ni, n2, and rir as the internal labels of their respective 
collecting edges. 

Although a rigourous derivation of the provenance of the angle operator is 
rather complicated, it is possible to gain an intuitive understanding as follows. 
The angle operator will measure the angle between the "bunches" of edges that 
traverse 5*1 and those that traverse 6*2. The most "natural" direction that is 
associated with the edges traversing ^i, 5*2, and Sr is the angular momentum 
vectors associated with their internal edges (which we define as Ji, J2, and Jr 
respectively). Classically, we would expect the angle between Ji and J2 to be 
given by 

cose^ (3.34) 



Ji 



J2 



^''Note that even if the internal structure of the intertwiner is not of this form for the parti- 
tion we choose, i t ca n be viewed as a superposition of such intertwiners (using the techniques 



stated in section 3.2.) 
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•J2 



-J3 




Sr 



C3 -s. 1 

Z3 ^ I ^ 



Z2 



Figure 9: Internal decomposition of a vertex operated on by the angle operator; 
the dotted hne denotes the boundary between external and internal edges. All 
of the edges xi^X2t ■ ■ ■ which intersect the surface Si arc "collected" into the 
internal edge labelled with n\; similarly, all the edges which intersect S'2 are 
collected into the edge labelled 112, and the edges which intersect Sr are collected 
into the edge labelled n^. In the text, we denote si = X^Xj, S2 = J^Vi' 
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^0 P, Pt Ps 



Figure 10: Internal decomposition of a vertex to be operated on by the volume 
operator. The external edges Pq, Pi, ■ ■ ■ , Pn-i are "lined up" along an internal 
"spine." The edges Z2,«3, . . . ,«n-2 are internal. 



However, since we don't want a net angular momentum to be associated with a 
point in empty space, we require that 

Ji + J2 + X = 0. (3.35) 

Using this, we can put the quantity in ( 3.34 ) in terms of operators: 

(X)2 = (Ji+J2)2 

= Jl + Jl + 2.h-J2 

Jl-J2 = l{Jr-Jl-J2) (3.36) 

The angle operator is then defined as 



e = arccos ( J^2 \ (3 



1\J J\ \jJ^ 



2 

Since we know the action of the operators already, we can readily state the 
eigenvalues of the angle operator: 

/ nr(nr + 2) - ni(ni + 2) - ^2(^2 + 2)\ 

d — arccos _ , (3.38) 

\ 2^ni(ni + 2)^712(^2 + 2) ) 

We note briefly that there are important philosophical considerations which 
arise from the discrete spectrum of this operator as it relates to diffeomorphism 
invariance. These considerations are discussed in Appendix we will not go 
into them here. 



3.5.4 The Volume Operator 

The volume operator also acts on a vertex. In this case, however, it is more 
convenient to use the basis pictured in Figure |l^, where all the edges are lined 
up sequentially. As we did for the area operator, we will also find it easier to 
define the volume-squared operator (known as first; the volume eigenvalues 
are then the square roots of the absolute values of the eigenvalues of W . 
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The action of W itself is rather comphcated. For a vertex with n edges 
labelled {Pq, Pi, ■ ■ ■ , Pn~i}, we define an operator W[rst] for each triad of edges 
such that 0<r<s<i<n — 1. The action of W^^st] consists of three 2-edges, 
joined at a trivalent intersection, which grasp edges r, s, and t as follows: 



(3.39) 



p, p, 



Pr P, 



= PrPsPt 



r 








i 

2 


2 





(3.40) 



Since the "comb" basis spans the space of all intertwiners, we can write the 
action of the W operator as sending the original vertex to a superposition of 
other vertices in the same basis: 



p, p, 



p„-i 



Po Pr P, 



E 

fc2,... 



rst]i 



(3.41) 



This is essentially a matrix equation; the entries of the matrix are now indexed 
by the internal intertwiner edges. We can find these entries through extensive 
use of the identities in Section |^. This result was first found in ; the matrix 
entries turn out to be 



[rst\i2---in-2 



-PrP.Pt 



k2 Pt fcs 
12 Pt h 
2 2 2 



^j22 cfc4 . . . .fc„_2 



Tet 


' Pr 


Pr Pq 

i2 2 _ 


Tet 


' Ps Ps 

_ h h 


ki 
2 




0(fc2,i2,2)0(fc3, 


l3,2)e{Po,Pr 


k2)d{k2 


,k3,Pt)0{k3,ki,P,) 



(3.42) 



This expression simplifies somewhat if we express it in terms of the Kauffman- 
Lins 6-j symbols (see Appendix ^: 



r{n) k2---kn-2 
' [rst]i2---in-2 



wn::r.r-^ = -PrPsPti *2 Pt ^3 >x'^fst!---s^ 



k-n-2 
in~2 





Pr 






Ps ki \ 


\ k2 


«2 




1 fcs 


13 2 / 



e{k2,k3,Pt) 



(3.43) 
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In both these formulas we have used the 9-j symbol, which is equivalent to the 
following spin diagram: 




^3 





r 






J 



(3.44) 



The eigenvalues of the volume operator are then proportional to the sum of 
the absolute values of the W-eigenvalues: 



V = 



\ 



E 



0<r<s<t<ri-l 



16 



W[rst] 



(3.45) 



The reader may have noticed that this definition for the W matrices will not 
work in the case where the valence of the vertex is three or four. Of course, 
the grasping action of the W operator will be the same in these cases; the only 
differences will be in the way we calculate the matrix entries. 

In the case of the trivalent vertex, there is only one possible grasping: 



W, 



rst] 



P0P1P2 2 



^0 n f 



(3.46) 



Since there is only one possible intertwiner for a trivalent vertex, the space 
in which we are operating is one-dimensional; hence, we will end up with the 
volume-squared eigenvalue times the original vertex: 



r P 

2[]' y = 1^(3) 



(3.47) 



This matrix entry evaluates to be 







Pi 




P0P1P2 < 


I? 


Pi 








2 





e{Po,Pi,P2) 



(3.48) 



However, the 9-j symbol is antisymmetric under exchange of rows; hence, if two 
rows are the same (as in this case), the 9-j symbol evaluates to 0. Thus, the 
volume associated with any trivalent vertex is zero (a result first found by Loll 

In the case of a 4-valent vertex, things are a little more complicated. As in 
the rt-valent vertex (for n > 5), we will end up with a superposition of basis 



34 



states: 



(4) 
[012] 



= P0P1P2 





" r 




r 

2 


I 

9 





(3.49) 



k 



The coefficients W^[oi2]i^ can be evaluated using similar metliods to tliose used 
to obtain ( |3.42 ); the result is 



[012] 







Pi 


"A 




P0P1P2 < 


1? 


Pi 










2 







6{Po,p^,mP2.P5,mhj.'i) 



(3.50) 



De Pietri |17[ also gives an explicit polynomial formula for the matrix entries 
for an arbitrary 4-valent vertex with incoming edges a, 5, c, and d: 



(4) t+e 
[012] t-e 



c(_l)(a+b+c+<i)/2 



32^t(t + 2) 
[(a + 6 + i + 3)(c + d + t + 3)(1 + a + & - i) 
(1 + c + rf - f)(l + a + t - 6)(1 + 6 + t - a) 
(l + c + <-rf)(l + d + f-c)]^. (3.51) 

where t + e and t — e correspond to the internal edge of the vertex, and e = ±1 
(in other words, the matrix entries are non-zero if and only if the two states 
differ by 2.) Since this expression is polynomial in form (i.e. doesn't contain 
any combinatoric expressions), it is much easier to work with. We will return 
to this formula in Section 5.3. 
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4 The Angle Operator 



Although the form of the angle and volume operators are known, many of their 
properties have not yet been examined. In this section and Section we describe 
the results of original research done into the properties of these operators. 



4.1 Spectrum of the Angle Operator 



Recall from Section 3.5.3 that the eigenvalues of the angle operator are given 
by the formula 

e = arccos { ^^ijr ^ - 3.{j. + ^ - 3.{j. + ^ (4,^) 
V 2Vji(ji + 1)j2(j2 + 1) ; 

We will find it more convenient to use the core edge colours ni, ^2, and 
instead of their associated spins; in terms of these labels, the angle becomes 

/ + 2) - ni(rii + 2) - ^2(^2 + 2) \ 

= arccos , = , (4.2) 

V 2v/ni(ni + 2)n2(n2 + 2) / 



where 



n 



Note that the quantities ni, n^, and are the intertwiner edges that "collect" 
the spins from each of the three surface patches; they are related to the total 
edge flux intersecting each patch by the relations 

4 < si, n2 < S2, and rir < Sr (4-3) 

where Si, S2, and Sr are the respective surface fluxes. 

The easiest way to visualize the angle spectrum is to plot its possible eigen- 
values versus the total core spin (defined as ut = Ur + ni + 71,2). Such a plot is 
shown in Figure ^ This plot shows two interesting features. First, we note the 
existence of fan-like structures, consisting of several sequences splitting off from 
a sort of "spine" . (Such a structure has been isolated in Figure [ij. Second, and 
perhaps more importantly, there seem to be very few eigenvalues corresponding 
to small angles. We will deal with the latter of these phenomena first. 



To obtain a small angle, we want the argument of the inverse cosine in (4.2) 
to be as close to 1 as possible: 

g ^ nrjur + 2) - ni{ni + 2) ~ 7^2(^2 -|- 2) 

2 v/^ir(mT2>II(^I7T2) 

However, the internal edge spins are subject to the triangle inequalities: 

2ni, 2n2, 2nr < ni + n2 + fir 
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Figure 11: Angle operator spectrum for increasing total vertex spin. 



To maximize Q, then, we wish to maximize its numerator and minimize its de- 
nominator while simultaneously respecting the triangle inequalities. This occurs 
when 

fT-T nx 
rir ~ — and ni — no ~ — (4.5) 
2 4 

where nx is as defined above At these values of ni, n2, and n^, the argument 
of the inverse cosine then reduces to 

Q-^^- (4-6) 

For any angle e, then, we have the relation 



e > arccos 



Note that while the value for rir is exact, the values of ni and 712 are only approximate; 
this is because ut must be even, but may not necessarily be divisible by 4, and because ni, 
n2, and rir must be integers. However, for large values of ny, the values of ni and n2 that 
maximize Q will approach ut/^- 
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Angle Eigenvalue Spectrum vs. Total Spin of Vertex 



90 



30 



10 - 
- 



1 1.5 2 

Angle Eigenvalues (radians) 



Figure 12: One of the fan-like structures that develop as m increases. This 
particular fan consists of sequences for which m + n2 — rir = 4. 
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For small e, cosx is a uniformly decreasing function; hence, for e <C 1, we have 



cose < 



< 1 



e > 



4 



VnT + 8 



(4.7) 



Hence, for large values of ut, the minimum possible angle between two 

— 1/2 

patches is roughly proportional to . Moreover, we note that if sr > 

si + S2, then tit < 2(si + S2); this implies that for a "macroscopic" vertex, 
where Sr 3> si, S2j the minimum observable angle depends only on the edge flux 
through the patches whose angular separation we are measuring. 

The first observation (i.e. the fan-like structures) is somewhat harder to 
quantify. A bit of investigation shows that the "feathers" of the fans in Figure 
|l2| are sequences of the form 

rir — m, ni = m ^ r , and 77-2 = g, 

where m varies and r and q are constant. From the triangle inequalities, we find 
that 

\r\ < q and r + q is even. 
The angles associated with such a vertex are given by 



cos^ 



2mr - 



2q 



2y/q{q + 2)(m2 - 2m(r - 1) + r(r - 2)) 

and for r,q<^ m, the sequences will approach the value 

r 

cos 6 = — , 



(4.8) 



(4.9) 



It is natural to ask whether any real number in the interval [—1, 1] is arbitrarily 
close to a number of this form; if this were true, then it would follow that any 
angle could be approximated by a vertex with sufficiently large internal spin. It 
turns out that this is true: 

Lemma 1 For any real number x d [—1, 1] cind any e > 0, there exist integers 
q and r such that 



< e. 



Proof: Choose integers b and c such that 



< 
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this can be done because the rational numbers are dense in the reals. Further, 
choose another integer n such that 



1 - 



v/1 + 2/nc 



< 



2' 



this can be done since the quantity on the left-hand side of the equation can be 
viewed as a function which approaches as n ^ oo. So, if we set r = nb and 
q = nc, we have 




< X 

< e 



Hence, we can approximate any angle 9 to arbitrary accuracy as long as we have 
a vertex of sufficiently high spin. We note that this is merely a confirmation of 
Moussouris' spin geometry theorem [p^ . 

It is important to note, however, that while this proof shows the existence of 
vertices which yield any angle to arbitrary accuracy, it does nothing to associate 
a minimum total spin to this angle and accuracy. As it turns out, this task 
would be fairly di fficu lt to accomplish by using these sequences. If we expand 
the expression in (4.8) to first order in m^^, we end up with 



COS0 = 



I 



2^/9(9 + 2) 
1 



(2rnr - - - 2g)(m^ - 2m(r - 1) + r{r - 2))"^/^ 

r - 1 



2m^q{q + 2) 



(2mr — r 



r-1 



2q)[l 



2q 



yJq{q + 2) \ m j 2my/q{q + 2) 

r ^ I r{r - 2) + q{q + 2) 
Vlil + 2) ^ 2^q[q + 2) 



(4.10) 



We see that in general the rapidity with which this sequence approaches its 
limit is very difficult to obtain: it depends non-trivially on the choice of r and 
g, which in turn depend non-trivially on the angle we wish to approximate and 
the accuracy to which we wish to approximate it. 

In the special case where x can actually be written in the form r/ \J q{q + 2), 
however, it is fairly simple to determine the minimum spin required: simply plug 
r and q into the second term of (4.10) and find the minimum m that reduces 
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this term below the desired accuracy. Given this fact, the natural question to 
ask is whether any real number x can be written in this form. Unfortunately, 
this is not the case: any number of the form r/ ^ q{q + 2) is algebraic, so if x 
is transcendental it cannot be written in this form. Let us examine a case that 
is even more specialized: suppose x is of the form ^JTifb^ where a and h are 
integers. We wish to find integer solutions to the equation 

o 



which is equivalent to the second-order Diophantine equation (in q and r) 

aq^ + 2aq - hr^ ^ Q (4.11) 

If we make the substitution p = q + \, this equation takes on the form 

ap^-br'^~a = Q (4.12) 

which has the (trivial) solution p = l,r = 0. As it happens, there are recursion 
relations for Diophantine equations of this sort jl^ ; if we know one solution to 
a given equation, we can find other solutions using these relations. In this case, 
if there exist integer solutions u and v to the equation 

- abv^ = 1, (4.13) 



then p = u,r = av is a solution to (4.12) 



■2 ba^v^ -a = a{u^ -abv^ ~l) = 0. (4.14) 



Equation ( 4.13|) is known as a Pell equation; it has solutions when ab is 



not a perfect square.]^ However, the solutions are not always small, even for 
fairly small values of a and b. For example, if a = 1 and b ~ 61, the smallest 
non-trivial solution to ( fl.ll| ) is: 

r ^ 226153980, q = 1766319050. 



Also, solutions to (4.13) are found by the expansion of v ab in a continued 

the forn 

ab = aQ 



fraction, i.e. an expression of the form 

1 



ai 



(12 + - ■■ 



Since the coefficients {ao, ai, a2, . . . } vary widely depending on the choice of a6, 
so do the solutions .|^ Hence, we are forced to conclude that finding any angle 6 
to an arbitrary accuracy e using this "fan" method is possible, but the solutions 
are difhcult to find and have an extremely complicated dependence on 9 and e. 

Note that when ab is a perfect square, y^a/b = a/ •/ab is rational. 
^^How these continued fractions relate to the solutions of the Pell equation is a fascinating 
subject which is, unfortunately, outside the scope of this piece. The interested reader is 
referred to Chrystal Eol. 
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4.2 Mean Angular Resolution 

Despite the difficulties in approximating a specific angle to arbitrary accuracy, 
it is possible to get an idea of the average separation between angles over the 
interval [0, tt]. There are two ways to look at the mean angular separation. The 
simplest way of finding the "angular resolution" for a given vertex is to note that 
an n-valent vertex has n—2 trivalent vertices in its intertwiner. Since each one of 
these internal vertices can be looked at as a "core" vertex, we can conclude that 
such a vertex has at most 3(n — 2) possible eigenvalues for the angle operator 
(the factor of three comes from regarding each of the three edges at each vertex 
as the rir edge.) We define the mean angular resolution 5 to be the average 
separation between each possible angle and the next greatest possible angle; in 
other words, 5 is the average width of a "gap" between angles in the interval 
[0, tt]. Since the 3n — 6 eigenvalues associated with the vertex divide the interval 
into 3n — 5 gaps, we conclude that the mean angular separation between the 
eigenvalues of the angle operator for an n-valent vertex is 

^ total of all gaps tt 
number of gaps 3n — 5 

However, this method does not take into account the fact that internal ver- 
tices may yield the same angles. For example, the internal vertices of a 4-valent 
vertex with all edges the same colour are identical; hence, when each of these 
internal vertices is treated as the core vertex by their respective angle opera- 
tors, there will only be 2 possible angles instead of the expected 6. Since this 
degeneracy depends on the internal structure of the vertex, it is impossible to 



say much about it in general; the best we can do is to modify (4. IE) to read 



This method, while legitimate, is somewhat limited: it only examines the 
separation between the eigenvalues of the angle operators that have a given 
vertex as their eigenstate. However, other angle measurements are possible for 
a given vertex. For example, if we measure the angle between two edges that 
do not immediately intersect in the vertex's intertwiner, this measurement will 
not always yield a definite eigenvalue; instead, there will be an expectation 
value associated with this measurement. The question then becomes: for a 
given vertex, can we find the total number of possible expectation values for all 
possible measurements? 

We note that given any vertex, we can partition its edges into three groups 
(one for Si, one for 5*2, and one for Sr-) If we switch basis to the "snowflake" 
basis, we will generally obtain a superposition of edges of this form; however, 
for a given superposition (corresponding to our original vertex) and a given 
partition, we will find a single expectation value. Hence, for a given vertex, the 
number of expectation values for the angle operator is equal to the number of 
partitions of its edges into three distinct groups. 
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Finding this number for a vertex with an arbitrary number of edges with 
arbitrary colours, however, is another question entirely. We note that the pos- 
sible range of colours for the collecting edge is not dependent on the order of 
the edges along the branch; hence, the number of distinct partitions is given by 

oo 

f^=Y[B{q,), (4.17) 

4=1 

where qi is the number of edges with colour i, and B{n) is the number of distinct 
ways of partitioning n into three distinct bins (one bin for each branch. 
Elementary combinatorics tells us that the number of ways to do this is 

5(n).r(" + ^)-^Ui!i±ll(!i±^ (4.18) 



However, we must divide the overall product by two, since exchange of the sets 
corresponding to the ni branch and the n2 branch yields the same angle. Thus, 
the maximum number of possible expectation values for a given vertex is 

i=l 

In the case where all the edges have the same colour, this formula reduces to|^ 

^^(n + l)(n + 2) 
4 

and the mean angular resolution (between expectation values) is given by 

47r 

S>^ (4.21) 



One might still wonder whether the expectation value of the operator could somehow 
depend on the order of the valent edges along the branches. To show that it does not, let us 
consider two states 

= ^ ai|ni,n2,n,-,Oi> and |'i/'2> = ^ fei |ni , "2, "r , ft> 

i i 

where ni,n2,nr are the collecting edges, and cxi and Pi are all other indices that can vary 
within the ni, 712, Ur subspace. (Note that the a indices and the /3 indices are not necessarily 
in the same basis, but both bases span the n\,n2,nr subspace.) The expectation values for 
each of these states are then 

(V'l|6'|t/'i> = /(ni,n2,n,.)^ |ai|2 and (V'2|e|V'2> = /("l, "2, "r) ^ 



where f{ni,n2,nr) is the function of ni,n2,nr given in (4.2). Since application of the recou- 
pling theorem is a unitary basis transformation, = '}2i l^iPi hence, the expectation 

values of both states are the same. 

^^In general, we will call a vertex whose valent edge colours are all the same a monochromatic 
vertex; we will also refer to a monochromatic vertex whose edges are a specific color as a 
monochromatic vertex with edge colour m, or sometimes as an m-chromatic vertex. 
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where we have again used an incquahty because of possible angle degeneracies. 

It is important not to confuse the "mean angular resolution" as we have 
defined it with the average distance from a point in the interval [0, tt] to the 
nearest angle (averaged over this interval). The quantity S that we have found 
allows us to get an idea of how tightly, on average, the eigenvalues of the angle 
operator for a given vertex are spaced. To find this other value (which, if 
found, would help us to estimate the level to which an arbitrary angle can be 
approximated) would 

f prob. of being \ /mean separation from\ 
V located in the gap / V either end of gap / 

all gaps ^ ^ ^ \ ^ , 

This is equivalent to the expression 



(^) = E 5 (4-23) 



47r 

a 

where Wa is the width of the gap a. This expression depends heavily on the 
precise location of the eigenvalues in the interval [0, tt] , and is therefore difficult 
to compute analytically. 



4.3 Angle Distribution 

In the classical continuum model of space, the distribution of solid angles is 
proportional to sm0: 

r{e)de = smede (4.24) 

If our quantized angle operator is to be of any physical worth, it must reproduce 
this distribution in some "classical limit." To find out whether this is so, we 
must examine not only the location of possible eigenvalues in the interval [0, tt] 
(as we have in the past two sections), but also the likelihood with which these 
angles occur. Due to the complexity of this problem, we will consider only the 
simplest possible case which could still conceivably have a classical limit: that 
of an n-valent monochromatic vertex with edge colour 1. 

We note that every vertex can be transformed into the snowflake basis by 
repeated application of the recoupling theorem. If we characterize these states 
by the internal spins a, b, and c, this fact can be written algebraically as 

IV') = ^Cabc\ai,a2,.-. ,61,62,- ■• ,ci,C2,...) (4.25) 

a,b,c 

where a, b, and c (without index) denote the sets {ai,a2, • ■ • , ai}, {61, 62, . . . , 6^}, 
and {ci, C2, . . . ,Ck} respectively. The most difficult part of this formula to an- 
alyze for arbitrary n is the coefficients Cabc! finding these for a specific internal 
intertwiner structure could be done, although such a calculation would be mind- 
bogglingly tedious. For a random vertex, however, we do not know the exact 
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internal structure a priori; hence, there is no reason to assume that any given in- 
tertwiner is preferred. We will therefore assume that all intcrtwiners are equally 
likely. With this simplification in mind, we can conclude that the probability of 
measuring the angle associated with an intcrtwiner core ni,rt2,«r is given by 

, , / number of intertwiners\ 
P0ni,n2,n, cx 4.26 
\ with core rii, 712, / 



where 6{ni,n2,nr) is given by ( |4.2| ). 

We wish to know, then, how many intertwiners exist with a given ni, n2, 
and rir- This is found in the following lemma: 

Lemma 2 If i is the number of edges ( of colour 1 ) entering a branch, the 
number of distinct branches with internal edges {02,03,... ,0^} that end in a 
given edge label ai is 

Proof: Consider the network depicted in Figure |l3|. Given a sequence of inter- 
nal edges {fli, 02, . . . , On}, we can define a path on this network which passes 
through the points {(l,ai), (2,02), . . . , (n,a„)}. Thus, the number of distinct 
branches with internal edges {ai, 02, . . . , a^} is equal to the number of distinct 
paths on this lattice from the point x — l,y — I to the point x — i,y = at 
(which go strictly to the right at all times, not to the left.) We also note that 
Q(i, Oi) — Q{i — 1, Oi — 1) + Q{i — 1, ai + 1); in other words, the number of paths 
to a given point in the network is equal to the sum of the numbers of paths 
to the two points closest to it in the previous column. Using this technique, 
the number of paths for small values of x and y can be easily calculated; these 
numbers are shown in Figure [13^, above each vertex. 
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"t 



J- = 3 



A- = 5 



Figure 14: The lattice in the m-n plane. To make the correspondence between 
this figure and Figure |l3| clearer, we have included two representative lines of 
constant x and y in this diagram. 



We now make the substitutions 

X^ry x-y 

m — and n — ; 

2 2 

note that m and n are integers since x = y (mod 2) for any point on our original 
lattice. If we redraw our lattice in the m-n plane, we obtain Figure Our 
paths are now on a conventional lattice; however, the condition that y > for 
the paths on our original lattice translates into the condition that m> n (since 
y = m — n.) Thus, we wish to find the number of lattice paths from (0, 0) to 
(too, no) (where too = (i + ai)/2 and no = (i — ai)/2) which never go above the 
line m = n. 

This new problem is a classic problem in combinatorics; it is discussed, along 
with other related problems, by Hilton and Pedersen [Q. The result for the 
lattice paths is then 

^1 . Too - no + 1 Atio + no + 1 

Q(TOo,no = ■ — - 

nio + no + 1 V nio + 1 

and the solution to our original problem is then 

^ a^ + lf i + l 
Q(t, ai) = 



1 + 1 



If we wish to turn this into a normalized probability distribution, we can 
approximate the binomial coefficient in the above expression by a Gaussian 
distribution multiplied by a normalization factor A: 



Aexp 



y\2-\ 
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Applying this to (4.27), we have 



P{h ai) = , . exp 
t + 1 

where A has been subsumed into the exponential term. 

Using this probability distribution, we numerically computed the distribu- 
tion of the angles in the interval [0, tt] . The program randomly selected candidate 
values of rii, n2, and Ur, making sure that they satisfied the vertex constraints; 
it then accepted or rejected the angle using the probability distribution derived 
above. The results for various values of si, S2, and Sr are shown in Figures ^ 
to H. 

Figure [T^ shows the distribution for si = S2 = s^. This distribution certainly 
does not match the expected classical distribution; most notable is the fact that 
the distribution peaks at approximately 2.5 « 140°. 

Figures |l^ and |l^ show the distribution for si — S2 Sr and si < S2 ^ Sr, 
respectively. These have surprisingly good correspondence with the expected 
classical angle distribution. 

Figure ^ shows the case Si <C S2 <C s^. This particular case required 
an extraordinary amount of computer time; as a result, the deviation of the 
data points is much larger. However, the data still matches the expected sin 9 
distribution very well. 

These data, while not conclusive, are certainly indicative that the angle 
operator can indeed reproduce the classical angle distribution in the case where 
si,S2 ^ 1 and si,S2 ^ Sr- This first condition corresponds to what normally 
thinks of as a "classical limit," well away from the regions where quantum 
effects dominate. The second condition can be thought of as a requirement for 
a sufficient amount of "background geometry" : not only do we need to make 
sure that the angles we are measuring include a large amount of spin, but we 
also need to ensure that the background spacetime upon which we measure the 
angles is sufficiently classical. 



2{i + l) 



(4.28) 
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Figure 15: Angle distribution for si ~ S2 — Sr- 




48 




Figure 17: Angle distribution for si < S2 <^ Sr- 
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Figure 18: Angle distribution for s\ <C S2 <C s^. 
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5 The Volume Operator 



5.1 The VF- matrix and Eigenvalue Bounds 



We recall from section |3.5.4| that the elements of the Ty-matrix are given by 

k2 Pt h 



^[rsf]i2-«„-2 - -^Ps^t 



*2 Pt H }^^s^:---si., 

2 2 2 



where fc2, 



, kn-2 and i2, ■ 



r Pr 


-P. 


Po\( Ps Ps ki \ 






2 / \ fcg »3 2 / 






e(fc2,A:3,Pt) 


5 *n2 


are 


the internal spins in 



(5.1) 



(see Figure ^.) 

From this operator, we then define the "volume-squared" operator as 



r<s<t 



[rst\ 



16 



(5.2) 



where is once again the Planck length. The eigenvalues of this operator are 
the squares of the volume operator V. 

We can see that this expression is horrendously difficult to handle in all its 
generality, for three main reasons. First, for an n-valent vertex with arbitrary 
edge colours, there are (3) possible choices for r, s, and i; each of these choices 

(n) 

could conceivably produce a different VFj^^^j. Second, the eigenvalues of the 

operator are difficult to find; even those of the W^-matrix are distinctly 
non-trivial, since all that ( |3.l[ ) gives us is the matrix entries. Third, the matrix 
entries themselves, as given by (5T), are not easy to calculate analytically. These 
three factors make the overall problem rather intractable, either analytically or 
numerically. 

With these problems, it is to our advantage to make some simplifications. In- 
stead of considering an arbitrary vertex, we will examine the case of a monochro- 
matic vertex; in this case, all of the matrices W^^^^,a,ie identical (since our 

"graspings" are always on edges of the same colour. )q The formula (5^) then 
reduces to 



G) 




n{n- l)(n-2) 




16 


96 



(5.3) 



This solves our first problem. The second problem, however, is somewhat 
less tractable; analytically finding the eigenvalues of an rt x n matrix is equivalent 
to finding the roots of a nth degree polynomial, for which no analytic solution 

^"We will often suppress the indices [rst], since they are always the same. In other words, 



when referring to the monochromatic case, we define VF'") = W, 



(") ^ 

rst] ' 
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exists for n > A. The matrix W^"'\ on the other hand, has a row for each 
possible intertwiner core - a number that is, in general, much greater than 4. 
To get around this, we will examine bounds on the eigenvalues instead of the 
eigenvalues. For any n x n matrix with entries a^ , any given eigenvalue A of 
this matrix satisfies 

n n 

|A| < max > \o,ij\ and |A| <max> |ai,| (5.4) 

In other words, the magnitude of every eigenvalue must less than the sum of 
the absolute values of some column and of some row.^ Using this, we can 
find an upper bound AI on the magnitudes of the eigenvalues of VF^") . In the 
monochromatic case, however, the eigenvalues Xwa of the I^(") operator are 
related to the eigenvalues Xy^a of the V-^ operator by 

n{n~ l)(n-2) 

'^V^a^ TT^ \Xwa\ (O-O) 

96 



as can be seen from (5.2). Hence, we can place a limit on the eigenvalues of the 
volume operator: 



Xva — V Ay2„ 



n(n-l)(n-2) 

— \^Wa\ 



96 



S f'-'^f'^' M ,5.6) 

The determination of M — the last of our three problems above — is, 
unfortunately, rather intractable. While explicit expressions for the 6-j and 9-j 



symbols in (5.1) exist, they are rather complicated. Simpler expressions do exist 
for the case of a 4-valent vertex; in this case, we can easily calculate the matrix 
elements explicitly. However, in general, we will have great difficulty doing so. 



The other two well-known results for eigenvalue bounds — namely, Gershgorin's theorem 
and Schur's inequality — turn out to be much less useful to us. Gershgorin's theorem states 
that for some i, 

|A - aiil < ^ \aij\ 

and a similar result for the rows of the matrix; ho weve r, the diagonal elements of the W^"^- 
matrix are all zero, which reduces this formula to ( ^.4[ ). Schur's inequality states that 

n n 
ct — 1 — l 

where are the eigenvalues of the matrix. However, it is much more difficult to obtain 
an analytic expression for the RHS of this equation, as it includes all of the elements of the 
matrix. For more information on these bounds, consult Bell bal . 
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5.2 Volume Eigenvalue Bounds for 4-valent Vertices 

In the case of a 4-valent vertex, De Pietri ||l^ has shown that the entries of 
the W^^"^ matrix for a 4-valent vertex with edge colours a, b, c, and d are given 
explicitly by the formula 

[{a + b + t + 3){c + d + t + 3){l + a + b-t) 
{I + c + d - t){l + a + t - b){l + b + t - a) 
{l + c + t-d){l + d + t-c)]K (5.7) 

where t + e and t — e correspond to the internal edge of the vertex. It can be 
shown that these elements are zero unless e = ±1; hence, each row or column 
will have at most two non-zero elements. To find M in this case, then, all we 



have to do is to find the maximum value of the polynomial on the RHS of (5.7) 
In the monochromatic case, a = b = c = d; we will call this edge colour m 
The polynomial in question then becomes 



Ti/(4) t+e 
[012] t-e 



(2m-|-i-K3)(2m-i+ l)(t+ 1)2 
32v/i(i + 2) 



If we think of this as a continuous function VF(i), we can show that this expres- 
sion is maximized at 



traa. = ^ ^^"^ ^ ^ + ^ V(2m + 2)^ - 16(2m + 2)2 + 16-1, (5.9) 

and the maximum value is, of course, Wit-max)- In the large-spin limit, where 
771 ^ 1, this maximum value scales as 

vn 

\Wraa.\ « ^ (5.10) 

The maximum absolute row or column sum (i.e. the quantity M) is at most 



twice this number; therefore, (5.6) tells us that for a 4-valent monochromatic 
vertex with edge colour m 3> 1, any eigenvalue Ay of the volume operator on 
this vertex will satisfy the inequality 

™3/2 

\y < ll^= (5.11) 

This is a very encouraging result when we consider it in relation to the area 
operator. The eigenvalues of the area operator are given by 
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Figure 19: The eigenvalue limit (shown as a solid line) obtained in (5.8) and 



(5.!;), compared with the maximum volume eigenvalue for 1 < m < 100 (shown 
as points). The lower dotted line is the result of a curve fit of the maximum 
eigenvalues to a curve of the form A = krr?!"^ ^ where k is the only free parameter 
of the fit. This fit yields a value oik^ 0.22155. 



where the values rrii are all the intersections of edges with the surface in question. 
For a 4-valent vertex with large edge colour m, we have 

\a oc m. (5.13) 



Together with (5.11), this implies that the volume scales no faster than J^l'^ ^ 
which is what we would expect from classical spatial geometry. 

We must ask ourselves, however, how "good" this bound is: do the eigen- 
values actually scale proportionally to m^^^, and if so, is the proportionality 
coefficient we have found (i.e. l/\/l2\/3) close to the actual proportionality co- 
efficient? Figure ^ shows the largest eigenvalue of the volu me o perato r (f ound 



by numerical calculation) compared to the bound stated in ( p.8[ ) and (5_^). As 
we can see, the bound fits the volume eigenvalues very tightly, diverging signif- 
icantly from them for n < 10 and becoming practically indistinguishable from 



the maximum eigenvalue for n > 50. We also show, in Figure 20, the ratio be- 
tween the eigenvalue bound we have found and the maximum eigenvalue. Note 
that this converges rapidly to 1. 
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1 10 100 

Edge Colour (m) 



Figure 20: The ratio of the eigenvalue hmit to the maximum volume eigenvalue, 
for 1 < TO < 200. 



5.3 Eigenvalue Bounds for ra-Valent Vertices 

We now turn our attention to the more general case of the n-valent monochro- 



matic vertex. In this case, we cannot use the simple formula in (3.51), but must 
instead use the more general form in (B.42) or (3.43). There are some simplifi- 
cations that can be made which aid us in computing these matrix entries and 
their bounds; however, a general analytic result is not forthcoming. Numerical 
methods, however, will prove fruitful. 

First, we recall that we are examining the special case where P r ^ P s = 
Pt — Pg = ni. We the note that the expression for the 9-j symbol in ( 3.42 ) and 
( |3.43 ) is fairly symmetric, and can be simplified. This symbol corresponds to 
the spin network 



k2 m fcs 
i2 m iz } = 
2 2 2 



' . 2 


r 


m 
m 

h \ J 




J 



(5.14) 



We can, however, move the graspings of the 2-edges around using equation 
(B.IO), and then close the resulting "triangles" using equation (B.ll). Doing 
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this, we obtain 



r k2 


m 
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1 , C 
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'i2 Tet 



k2 12 «2 

2 2 2 



^(fc2,i2,2) 
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2 2 2 




«3,«3,2) 



Tet 



k2 fcg 2 
is *2 w 



Plugging this back in, (|3.42| ) becomes 



in) k2---kn-2 



[rst]i2---'in-2 



Tet 
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Tet 


ni 


m 




Tet 


'k2 


ks 
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k2 






ks 
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2 




. *3 


«2 


TO 



0(^2, «2, '2)9{k3, 13, 2)0{m, to, fc2)0(fc2, fcs, m)0{k3,k4, to) 



where 



TO 



-^2 Tet 



«2 «2 

2 2 2 



V 



e(fc2,Z2,2) 



igTet 


^3 h h 

2 2 2 




«3,«3,2) 



(5.15) 



We can also put the Tets in this resulting equation in terms of Wigner 6-j 
symbols (see Appendix |^) ; at least one of the arguments of the resulting Wigner 
6-j symbols will be 1, and analytic formulas for this case have been compiled 
(see, for example, Varshalovich [^.) Equation ( 3.42| ) then becomes 



W, 



(n) k2---kn-2 
rst]i2---'in-2 



(5.16) 



where Q is given in ( |5.15| ), and TZ and S are given by 



to/2 fc2/2 fc3/2 

1 «3/2 12/2 



w 



to/2 to/2 A:2/2 
1 12/2 to/2 



ki/2 to/2 fc3/2 
1 13/2 to/2 



w 

(5.17) 



<^_n( 9^ / 6>(^2,^3,TO)6>(TO,m,^2)6'(^3,fc4,TO) 

Finally, we note that there are certain "selection rules" for finding the non- 
zero elements of W in this basis. The presence of a series of Kronecker deltas 
in ( 3.42| ) gives us the obvious selection rules that 

^4=14, fc5=i5, kn-2^in-2 (5.19) 
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It is also evident from the spin diagram in (3.44) that for the 9-j symbol in 
( ^.42 ) to be non-zero, we must have 



A:2-i2 = 0,±2 and fc3-i3 = 0,±2 (5.20) 
in order to satisfy the triangle inequalities. Finally, it can be shown that 

2 



0(a,6,2) 



Tet 



a a b 
2 2 2 



(_l)«+i.ii±^ b 



-1 




a 

b — a 
b = a + 2 
otherwise 



(5.21) 



Hence, if k2 = i2 and fcs = 13, the quantity in parentheses in (5.15) will be 
(--1 + 1), with the result that Q = 0. We can then rewrite (5.20) as 



fc2-i2 = 0,±2, A:3-i3 = 0,±2, and 
fc2 — ^2 & ^3 — *3 are not both zero. 



(5.22) 



We conclude that for a given column of the W matrix (i.e. fixed fc2,/c3, 
fcn-2), there will be at most eight non-zero elements: 



«2 = ^2 - 2 
«3 = fcs - 2 

«2 = k2 



«2 
^3 



fc2 
fc3 



»2 = ^2 - 2 
«3 = ^3 



i2 = k2 + 2 



»2 = ^2 - 2 
«3 = fc3 + 2 

«2 = k2 

is = fc3 + 2 



(5.23) 



i2 

is 



k2 



Unfortunately, this is about as far as analysis will easily take us. The quan- 
tities Q, TZ, and S are fairly simple to calculate for a given set of arguments; 
however, finding the maximum product of all three, over all possible values of 
{^2, A:3, fc4, i2, *3j ^4} is a much more daunting task.Q Things are complicated 
by the fact that both the quantities Q and TZ have slightly different analytic 
forms for each of the eight entries in a given column — one form for fc2 ~ *2 = 2 
and fc3 — 13 = 2, one form for /e2 — Z2 = 2 and ks — is — 0, and so forth (see 
Varshalovich ^^.) We see, then, that an analytic answer is difficult to find. 

Despite this difficulty, we can get an idea of the scaling properties of the 
absolute row and column sums by numerical analysis. By a happy coincidence, 
we do not even need to look at all of the rows and columns of the matrix — 
only a small subset of them. This is because (as mentioned above) the entries 
of the ly-matrix are only dependent on the six values {fc2, A:3, fc4, Z2, 43, 14}; all 
other pairs of corresponding indices must be equal, or the matrix entry will be 



Note that, in a small mercy granted to us by equation (5.1), the entries of this form of 
the W^-matrix depend only on these six intcrtwiner strands (and the edge colour m), and not 
on the entire set. 
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zero. These facts imply that the W^-matrix must be of block-diagonal form, and 
that (more importantly) each one of the blocks is identical: 













M 
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w = 








A 







u 





• 


• a) 



(5.24) 



where each block corresponds to a different choice of k^, . . . ,kn-2, and the 
indices of the matrix A are all possible choices for fc2, k^, and /c4. We see that 
we need only compute the absolute row and column sums for the (much smaller) 
matrix A, instead of the full VF-matrix. This is still a non-trivial calculation, 
as the size of the matrix A scales fairly quickly (for m = 10, the matrix has 

already grown to 891 x 891,) but it does simphfy the pro blem greatly. 

Using this trick, combined with the selection rule s in ( |5.19| ) and ( 5.22 ), the 
bounding quantity M (as defined above in section |5.2| ) was calculated with 
MATLAB. Two methods were used. The first was the absolute row and column 
sum discussed previously. The second was MATLAB 's internal normest func- 
tion, which estimates the largest eigenvalue of a sparse matrix.^ These results 
are shown in Figure Each of these data sets is fitted to a power law rela- 
tion. We see that for both of these data sets, we have an approximate relation 
M oc m^. 

There is good reason to believe that this power dependence should tend 
towards (exactly) as m ^ oo: if it docs, then we have the relation 



|Av| 



n(n-l)(n-2) 



M 



96 



= C(total spin) 



3/2 



(5.25) 



where C is a proportionality constant. As noted in the previous section, this 
is the dependence expected from classical spatial geometry. While our data 
are not conclusive evidence of this relation, they are highly suggestive that this 
relation holds for a monochromatic vertex of any valence. 

Thus, we have placed limits on the volume eigenvalues of monochromatic 
vertices; in both cases, we find that volume of a region scales as the 3/2 power 
of the area of its bounding surface. This result was found analytically in the 
4-valent case; numerical results strongly suggest that this relation holds for the 
general case as well. 



For the mathematically inclined reader, we note that the absolute row and column sums 
are known as the 1-norm and the oo-norm of a given matrix, respectively; that which normest 
finds is the 2-norm. 
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^ — "Normest" matrix norm 
^ — Abs. row/col. sum 




1 10 

Edge colour 



Figure 21: The eigenvalue limits found by the absolute row and column sums 
and by MATLAB's internal normest function. The curve fits are for a power 
law relation. For the absolute row/column sum method, the best fit is M = 
0.90432 • m2-9i26. ^j. ^j^g normest method, it is M = 0.71884 • m^-^^^i^ 
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6 Discussion 



We have seen some interesting consequences stemming from the discreteness of 
the angle operator (see S ectio n 0) ; these resuhs wiU be referred to as the sma ll 



angle property (Equation 4.7), the angular resolution property (Equation 4.21 ), 
and the angular distribution property (Figures [T^-ll^.) We now wish to assign 
some scaUng behaviour to these predicted phenomena; in other words, on what 
length scales do we expect these phenomena to be observable? 

The small angle property is the best-quantified of these three properties; we 
have the equation 

4 4 
£> , (6.1) 

^/nT + 8 y/ST + 8 

where St is the total surface flux. From our numerical studies of the volume 
operator, we also have the result that 

V oce^ (X s^^ (6.2) 

where £ is the length scale (as defined by either the volume or area operator — 
we have shown that these are equivalent.) Hence, if we combine these results, 
we obtain the result that for large st, the minimum observable angle is related 
very simply to the length scale observed: 

e oc oc Srp^^^ (6-3) 



Moreover, from our discussion in Section 4.1, we recall that for a vertex with 
large amounts of "background geometry" (i.e. S> 81,82), the primary limit 
on ut is not st but si + 82- Hence, for these vertices we have the result that 

eoc (si +S2)"'/' (6.4) 

Using these relations, we can also gain an idea of the absolute length scale 
on which this phenomenon might be viewed. Recalling that 

y = l'o^/T^' (6.5) 

where M was found to be equal to 0.71884 (numerically), we can state an 
absolute relation between the length scale and the smallest observable angle: 

e = kjiy (6.6) 



where fee = 4 ^ AI/96 « 1.7692. This implies that if we observe angles on 
a length scale of ^ « 10^^^ m (i.e. the radius of a proton), the minimum 
observable angle is then approximately 10"^*^ radians — a miniscule angle by 
any account. 
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In Section 4.2, we found results for the mean angular resolution associated 
with a vertex. The most important of these was the average spacing between 
expectation values, given by 

477 

(6.7) 



Again, combining these results with the result in (6.2) allows us to associate a 
scaling behaviour with the mean angular separation and resolution. For large 
n, the angular separation scales as 

6(xn-^(xe-*, (6.8) 

This scaling property of the angular resolution is also somewhat striking, since 
it scales proportionally to such a high value of the total spin (i.e. inversely 
proportional to the fourth power of the length.) 

Using similar techniques to those found above, we can find absolute length 
scales for 6 as well; these can be shown to be 

S = ks (f) (6.9) 



where ks — 47r(M/96)^/'^ ~ 0.4810. On the length scale of a proton, this implies 
a resolution of roughly lO^^*' radians! 

Finally, we found that we could reproduce a classical angle distribution under 
the condition that 



l<{si,S2}<Sr. (6.10) 

This is somewhat harder to quantify; however, we can associate lengths with the 
distributions shown in Figures p^|-p^. For Figure we have tit ~ 10®; hence, 
the length scale associated with such a vertex is on the order of 10^ « 10^'^^ 
m; for Figures p^-p^, the length scale would be on the order of 10^° ~ 10~^^ 
m. For reference, the diameter of an atomic nucleus is on the order of 10~^^ 
m, a full 17 orders of magnitude larger than the length scale associated with 
the "drastically altered" distribution of Figure |l^. Even on a length scale ten 
orders of magnitude less than the size of a nucleus (as in Figures p^|-[T^), the 
classical angle spectrum is reproduced. 

We must be careful, however, that we do not read too much into these dia- 
grams. The main difference between the conditions for the normal and skewed 
angle distributions is that si,S2 ^ Sr for the former while si,S2 ~ s,. for the 
latter. To gain a better understanding of where the transition between these 
two types of distribution actually occurs, we would have to examine the distri- 
butions over a wide range of values of Si/s,. and S2/ Sr', the results presented 
here only show two values of si/sr and three of 82/ Sr (not counting the minor 
difference between Figures |l|and|l^). 

Moreover, while the skewed angle distribution in Figure |lj certainly occurs 
for si, S2, ~ 10^, there is no reason to believe that this behaviour will go away 
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for larger values of the spins. If this behaviour does indeed persist at higher 
spins, our interpretation of the angle operator may be incorrect. We recall that 
the angle operator finds the angle between the two patches Si and ^2; perhaps 
all this distribution is tolling us is that these patches don't "want" to overlap; 
as they get larger, they will have to spread themselves out to avoid each other. 
However, this is purely conjectural; more research is needed to find out whether 
this problem even arises. 

The true importance of our angle distribution results becomes clearer when 
we consider them in light of the small-angle phenomenon. We recall that the 
smallest angle observable about a vertex depends on st if si,.S2 ~ ^r, but on 
Si + S2 if si,S2 Sr- To obtain a classical angle distribution, the latter of these 
conditions must hold; the magnitude of Si and S2 then determines the minimum 
observable angle. Thus, if we have a vertex of a given spin, we have to balance 
two factors to obtain a classical angle distribution: the need to assign enough 
spin to Sr to obtain the correct distribution, along with the need to assign 
enough spin to si and S2- Hence, to approximate classical angles in both of 
these ways, the total spin of the vertex must be sufficiently high. 
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7 Conclusion 



The search for a quantum theory of gravity is one of the major chaUenges 
facing modern physics today. One of the more promising approaches in this 
search is the canonical or loop quantum gravity approach; in this theory, the 
fundamental variables of the theory are not field variables, but path integrals 
of field variables around closed loops in space. A difficulty arises from the 
properties of these loops: a given set of loops may satisfy any number of linear 
relations between themselves, making it difficult to compare two arbitrary states. 
This difficulty can be resolved through the use of spin networks, whose properties 
were described in Section 

Since general relativity predicts the structure of space, these spin networks 
must do so as well. This is done through Hermitian operators on spin network 
states, corresponding to "spatial observables" such as area, volume, and angle. 
While the forms of these last two operators are known, their properties have 
not been extensively studied. 

Our research into the properties of these operators has revealed several pe- 
culiarities. In particular, we have shown that the minimum observable angle 

— 1/2 

about a vertex with total spin st scales roughly clS 6 OC Srp , cl fairly slow de- 
crease to zero. Moreover, since the angle operator has a discrete spectrum, the 
possible observable angles will be separated from each other. We have shown 
that the average separation between adjacent angles in the spectrum decreases 
as (5 OC s^'^ . Finally, we have shown that for "sufficiently classical" vertices 
(i.e. I ^ {si, S2} ^ Sr), the classical angle distribution of 1^(9) d9 — sin^d^ is 
reproduced. However, in the case where si,S2 ~ Sr, the angle distribution is 
vastly skewed from its normal behaviour. 

Our studies of the volume operator were hindered by its complexity. How- 
ever, we have found (through a combination of analytical and numerical meth- 
ods) evidence suggesting that the volume associated with a vertex is propor- 
tional to the 3/2 power of a surface surrounding that vertex. This is identical 
to the scaling properties of these quantities in the classical model of space; 
hence, this result is a major step towards showing that spin networks can fully 
approximate classical continuum space. 

The true measure of any theory, of course, is not how well it predicts already- 
observed effects but whether it predicts any new, hitherto unexpected phenom- 
ena. In this sense, our angle results are more important than our volume result. 
The latter result is still important, of course; if we had found any other de- 
pendence than (volume) cx (area)'^/^, we would have had to discard our current 
volume or area operator — if not the entire theory — as untenable. However, 
the fact that our angle results predict unexpected phenomena (the small-angle 
phenomenon and the skewed angle distribution) makes these results rather more 
interesting. 

These phenomena imply that the angle operator predicts a drastically dif- 
ferent structure of space for low-valence vertices. Unfortunately, such vertices 
are very far removed from our current realm of observation, even if we consider 
the next generations of observational equipment. To get a rough idea of how 
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far away we are from experimentally probing the predicted properties of spin 
networks, we can note that to do so would require a particle with a de Broglie 
wavelength on the order of the Planck length. Hence, such a particle would have 
a momentum 

po ^ « 3 X 10^* cV/c « 16 kg • m/s. (7.1) 
to 

While this may not seem like terribly much momentum on a macroscopic scale, 
it is a tremendous amount of energy for an elementary particle to have. If 
we make the assumption that the rest mass of such a particle is very small 
(i.e. mo ^ po/c ~ 50 /ig), its energy will be (to first order) determined by its 
momentum only: 

= Poc « 3 X 10^* cV w 5 GJ. (7.2) 

This is a stupendous amount of energy; for comparison, the kinetic energy of a 
cruising Boeing 747 is roughly 12 gigajoules. We cannot possibly hope, then, to 
probe these distances directly at any point in the near future. 

Because of this, the skewed angle distribution result is rather interesting: 
it predicts a change in angle distributions at a scale three orders of magnitude 
larger than the Planck length. This would bring down the required energy of 
a particle required to see this phenomenon to approximately 10~^^ eV — still 
well beyond our observational capability, but a good deal closer. However, as 
mentioned in Section ^, this result is still open to interpretation. The normal 
distributions shown in Figures p^|-p^, however, do place an upper bound on 
the length scales on which quantum angle effects can be observed — by the 
time we have reached a length scale of 10^^'°^ m, the quantum angle spectrum 
approximates the classical angle spectrum very well. 

In the face of this, the small-angle property may end up being the "best" 
prediction stemming from our results. This is not to denigrate the other results 
we have found; the angle distribution and volume results show that the angle 
and volume operators do satisfy "classical" relations on a large scale. These are 
important properties for an quantum theory of general relativity to have, and are 
important steps towards showing that loop quantum gravity is a self-consistent 
example of such a theory. 
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Appendices 



A Dirac-Bergmann Constraint Analysis 
A.l The Generalized Hamiltonian 

To canonically quantize a classical theory, we need first to find a general ex- 
pression for its Hamiltonian. It is fairly simple to find one Hamiltonian for a 
given system; however, our objective is to find an expression for all the possi- 
ble Hamiltonian of a given theory. Such an expression should permit arbitrary 
variations that are implicit in the theory (e.g. gauge freedom in classical elec- 
trodynamicsO The following derivation of the generalized Hamiltonian follows 
from Dirac ||25|l . 

To construct the generalized Hamiltonian, it is easiest to start from the point 
of the action integral: 

5" = J L{qi,q2,... ,qN,qi,q2,--- ,qN)dt (A.l) 

where {qi,q2,... ,qN} are the coordinates of the configuration space of the 
system, and {gi, (72, ■ ■ ■ , qN} are their time derivatives. By applying the calculus 
of variations to this integral, we obtain the Lagrangian equations of motion: 

dqn dt ( dqn ) *^ ^ ^ 

Note that there are N such equations, one for each coordinate of the configura- 
tion space. 

The Hamiltonian formalism first requires us to change variables from the 
velocities g„ to the conjugate momenta p„: 

Pn = 1^ (A.3) 

dqn 

These momenta will satisfy the Poisson bracket relation 

{gn,Pm} = <5nm, (A. 4) 

where the Poisson bracket of two quantities / and g is defined as 

if . df dg df dg 

dqn OPn OPn Oqn 

It is possible that these momenta are all independent of one another; however, 
this is not required in general. As a simple example, consider a disc rolling 
without slipping on a plane. The disc's linear momentum px is not independent 
of its angular momentum pg . 

There will therefore be a certain number of constraints on the coordinates 
and momenta. These constraints can be written in the form 

(j)„,{q,p)^0 m = l,2,...,M (A.6) 
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where q and p denote any combination of the coordinates and the conjugate 



momenta. The expressions in (A.6) are called the primary constraints of the 
formalism. 

Let us now consider the quantity H = PnQn — L. We first note that an 
arbitrary variation in this quantity can be expressed as 

SH = Spnq, 



PnOQn - t; dQn - TT^Oqn 

Oqn Oqn 



A ■ ^^A 

oqn 



(A.7) 



Thus, a variation in H can be expressed solely in terms of variations in p„ and 
Qn] this implies that H is solely a function of the coordinates and their conjugate 
momenta, and not of their velocities qn- 

This quantity H is called the Hamiltonian of the system. However, it is 
not uniquely determined. Since the primary constraints are identically zero, 
we can add any combination of them to the Hamiltonian without changing the 
Hamiltonian's value: 



H* 



(A.; 



where the coefficients Um are arbitrary. (The negative sign is chosen merely for 
notational convenience later.) Using this new Hamiltonian, we can then deduce 



the generalized Hamiltonian equations of motion from (A.7) and (A 



dH* 

dpn 

dH* 



Qn 



dL 
dqn 



dpn 



du„ 

dpn 



Noting that the primary constraints vanish and applying (A.2) to the second 
equation, these equations lead to the Hamiltonian equations of motion: 



dH 

dpn 

dH 



dpn 



dqn 



(A.9) 
(A.IO) 



We can rewrite these equations in a fairly simple form using Poisson bracket 
notation. The time evolution of any quantity g{q,p) is given by 



9 



dg . 

^ — qn 



dg . 

dpn 



(A.ll) 
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Applying (A. 9) and (A. 10), we see that 
9 - 



dg /dH 



dg / 


'dH 

— -\- Urn 
^dpn 


d(t>ni\ 


9g„ ' 


dpn ^ 


dg 


( dcj)m , 




dqn 


^ dpn 


dpn 



dg dH dg dH 

dq„ dpn dpn dqn 

, dg ( d(j)jy, 

dPn ^ dqn dqn 

dg d{H + Ufa (pm 

) dg d{H + U,n(t>m) 



dur. 



dqn dp„ 
{9,Ht} 



dpn 



dqn 



(A.12) 



where Ht = H + Um(t>m is the total Hamiltonian. Note that the extra terms 
introduced in the second step are identically zero. 

This formulation of the Hamiltonian equations of motion is concise; however, 
it does pose a potential pitfall. Our primary constraints (jjm are, by definition, 
equal to zero; however, if we were to apply this equality immediately, before 



working out the Poisson brackets in (A.12), we would lose the generality that 
they allow. To clarify this distinction, we introduce the concept of "weak equal- 
ity", denoted 







(A.13) 



This notation signifies that one may only make use of this relation after the 
Poisson brackets have been worked out. Using this notation, then, ( A.12| ) be- 
comes 



{9,Ht}. 



(A.14) 



We can u se (A.14 ) to check for further consistency conditions of the theory. 
If we apply (A. 14) to one of the primary constraints 0™, we have 



'/>m ~ « {<j),n, H} + Um'{(t>m, 4>7n'} 



(A.15) 



Thus, for each value of m, we have a consistency condition. Each one of these 
conditions leads to one of three results :ri 



• The equations in ( A.15 ) could lead to an equation that is trivially satisfied, 
i.e. of the form = 0. 

• The consistency conditions could reduce to an equation that contains only 
coordinates and momenta. Such a condition must be independent of the 
primary constraints (otherwise, it would be identically zero and therefore 
be of the first type.) We therefore have a new set of secondary constraints: 



(l>k{q,p) ~ 



fc = 7W + l,X-|-2, 



(A.16) 



^■^It is possible that one or more of tliese equations lead directly to an inconsistency (i.e. an 
equation of the form 1 = 0); in this case, our original Lagrangian equations of motion were 
inconsistent, and the theory must be modified or rejected. 
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• Finally, the consistency conditions may not reduce in either of the above 
ways; in this case, they impose a constraint on the unknown functions Um- 

Secondary constraints may themselves generate other constraints by further 



apphcations of (A. 14), i.e. 

{<j)k,H} + u„^{<j>k,(l)m} ^0, (A.17) 

where k > Ai. These new constraints may be of any one of the three types 
listed above. If they turn out to be new secondary constraints, we can continue 



to apply ( A. 14 ) recursively to them. In the end, we are left with our original 
set of primary constraints 0™, a new set of secondary constraints 0^, and a set 
of equations of the third type relating the u„'s. 

Let us further examine these equations of the third type. Each of these 
equations is of the form 

{0j,i/} + u™{0j,0™}«O j = l,2,...,/C (A.18) 



We can consider ( A.lq) to be a system of /C linear non- homogeneous equations. 



with the Urn's as variables and coefficients in terms of the commutators of the 
primary and secondary constraints. Since these coefficients are in terms of p's 
and q's, we conclude that the u^'s must also be functions of the coordinates 
and the momenta: 

Um^Umip,q) (A. 19) 

A complete solution for these variables also includes an arbitrary linear com- 
bination of the independent solutions (if any) to the homogeneous system of 
equations 

K«{</'j,</'„J = 0. (A.20) 
Hence, the general solution to this system of equations is 

U„i = U,n + VaVarn, (A. 21) 



where a sums over the solutions found in ( A.20 ). The total Hamiltonian can 
then be written as 



where 



and 



Ht = H'+ Va(pa (A.22) 
H' = H+ Um(l>m (A.23) 
4>a = Vam4>m- (A. 24) 
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It is important to note that while the original Um coefficients were not neces- 



sarily independent of one another, the new coefficients Vm in (A.22) are totally 
independent .p^ 

At this point, we introduce a bit of nomenclature. We call a dynamical 
variable R first-class if its Poisson brackets with every constraint (primary and 
secondary) are weakly zero: 

{i?,0,}«O (A.25) 

An expression that is weakly zero must be equal to some linear combination of 
the constraints; hence, we have 

{i?,0,} = c,y0,. (A.26) 



If R does not satisfy ( A.25 ) for each </ij, we say that R is second-class. The 
"class" of a constraint is independent of whether it is a primary or secondary 
constraint; hence, we now have a total of four types of constraints. 



Let us again consider the arbitrary coefficients Va in (A.22). These new 
coefficients enter into the equations of motion: 

g-{g,H' + Va<l>a} (A.27) 

However, a classical theory such as ours should be deterministic; we must there- 
fore inquire as to the meaning of these arbitrary coefficients in the equations of 
motion. 

Initially, the state of the system is determined solely by the coordinates 
and momenta, not by the VaS. If the theory is deterministic, then this state 
uniquely determines the system's later state. However, since the evolution of the 
coordinates and momenta depends on the values of the coefficients Va, certain 
values of coordinates and momenta must correspond to the same physical state. 
To determine the relations between these values, we consider an infinitesimal 
time evolution in any dynamical variable g: 

Stg = .90 + gSt 

= go + {g,HT}St 

= go + St{{g, H' } + Va{g, (A.28) 

The Va coefficients are completely arbitrary; thus, we change them to obtain a 
new Stg. Since our theory is deterministic, these two new ^'s must correspond 
to the same state; therefore, the difference between them must correspond to a 
transformation that does not change the state: 

Ag = St{va - v'a){g,(t)a} 

= ea{.9,0a} (A.29) 



'In general, there will be fewer of the Vm's than of the Um's. 
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We can therefore change g using this rule, and the physical state will remain 
the same. These are known as infinitesimal canonical transformations^ and the 
0a 's are their generators.^ 

Suppose we apply two such transformations successively. We obtain 

9 ^ 9 + ea{9,(l>a} 

^ 5 + ea{g,0a} + 7Q'{.9 + ea{5,0a},0Q'} = 5' (A.31) 

If we apply them in the opposite order instead, we get 

.9'~^.9 + 7a'{5,0a'} + ea{.9 + 7a'{ff,0a'},0a} =9" (A.32) 

Since these two values of 9 correspond to the same physical state, a transfor- 
mation corresponding to their difference must also be an infinitesimal canonical 
transformation that is not associated with a change of state: 

9' - 9" = ea{5,0a} +7a'{5,0a'} +7a'ea{{g,0a},0a'} 

- 7a'{5,0a'} - ta{9Aa} " l^ala' {{9 Aa'} Aa} 
= eQ7a'({{5, 0a}, 0Q'} + {{0a', g}, 0a}) 

= W ({5, {0a, 0a'}}) (A.33) 

We see that the commutators of the 0a 's can also generate infinitesimal 
canonical transformations. However, it is important to note that while the orig- 
inal set of canonical transformations were generated by primary constraints, the 
new transformations generated by the commutators are not necessarily primary 
constraints. 

This leads us to the extended Hamiltonian He, which consists of the total 
Hamiltonian Ht plus any linear combination of expressions 0a' which generate 
infinitesimal canonical transformations not associated with a change of state: 

HE=HT+V',4a' (A.34) 

We will find in the next section that in the case of electrodynamics, these canon- 
ical transformations will correspond to gauge transformations. This last equa- 
tion, then, is the most general expression for the Hamiltonian of a system. 



A. 2 Dirac-Bergmann Analysis &z Electrodynamics 

As an example of Dirac-Bergmann analysis, let us attempt to construct the 
generalized Hamiltonian for classical electrodynamics. We saw in Section 
that the field tensor is given by 

F^, ^ d^,A, ~ d,A^. (A.35) 

^^For our purposes, we will consider an infinitesimal canonical transformation to be a trans- 
formation that adds an arbitrary infinitesimal (possibly a function of the original variable) to 
9, e.g. 

g g + ta{g,(pa} (A.30) 

In the above equation, (pa is called the generator of the infinitesimal contact transformation. 
For further explanation of canonical transformations, consult Goldstein |26|. 
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In terms of this tensor, the action can be shown to be 

S[A^] = -\jF^'F^,A^^ (A.36) 

where the integral is taken over a region of spacetime. 

Constructing the momenta from the action requires a Lagrangian; normally 
this Lagrangian is related to the action by the relation S = J Ldt. To obtain 
our Lagrangian, then, we must define a reference frame; this defines a f-direction 
on spacetime. We can then split the four-dimensional integral in (A.36) into a 
double integral over space and time: 



5 = -- // F^^F^.d^xdt (A.37) 



1 

4 

The Lagrangian is then the spatial part of this equation: 

L = -'^jF^-F^,d?^ (A.38) 

In a simple particle-dynamics system, we find the momenta by taking the deriva- 
tive of the Lagrangian with respect to velocities, as in ( |A.3| ). Taking the deriva- 
tive with respect to a field, however, is less straightforward. 

To find these momenta, we first define the Lagrangian density as a scalar 
field over all space: 

£(x) = -V''(x)F^.(x) (A.39) 

Note that this is simply the integrand of ( A.38| ). The momenta, which we rather 
suggestively label E^^ ^ are then defined by 

E^{^) = ^^i^ (A.40) 

To find these explicitly, we split the Lagrangian density into time-dependent 
and non-time-dependent parts: 

C = -i(^^™Foo + F'^'^Fao + FO-Fo, + F'^^Fab) (A.41) 

The first term of this expression vanishes, since i^^^ is anti-symmetric. More- 
over, the last term does not contain any terms containing the velocities (i.e. 
terms of the form d^A^j^), so we can discard it when taking the derivatives: 



dA^ 



= -l^iF'^'Fao + F^^Foa) 
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Recalling that A^^ = doA/^ and that = —F°-^, this reduces to 

Ei^ = F^o. (A.42) 



Wc sec now why the symbol E'' was used to denote the momenta: the 
momenta turn out to be the components of the electric field. Moreover, since 

= F^^ = 0, we can discard the t-component of this vector; hence, we will 
also denote the electric field as E"-. Moreover, since F" = 0, its conjugate 
variable Aq is not physically significant (i.e. Aq is a cyclic variable); hence, we 
will denote the configuration variable well. However, it is important to 

note that these constraints are only weakly equal to 0; to retain all generality, 
we will have to retain E° and Aq until we have found the secondary constraints. 

The Poisson brackets of the configuration variables Aa and E" will then 
satisfy the relations 

{E'^ix),E''iy)}=0 

{Aaix),Abiy)} = (A.43) 
{A4x),E\y)}=5lS\x~y) 

Now that we have the momenta, we can construct the Hamiltonian for the 
electromagnetic field: 

H = j E^'^oA|, d^x - L 

= j F^O^o^a + \ (i^°°F00 + i^"°i^aO + i^°"^ba + F^'^Fab) d^x) 

= j F'^^daAo - F„o) + iF«°F„o + ^F-'Fab d^x 

where wc have used the relation doAa = daAo — Fao in the first term and 
collected two terms together to make the second (through the symmetry of F.) 
Rearranging and combining terms, we have 

H = j ^F'^'Fab - ^F-^OFao + F'^^'daAo d^x 

= J iF«^F„6 + ^F«F„ - AodaE- d^x + J AoE^Ua d^x (A.44) 

Note that we have applied an integration by parts in this last step, creating 

the last term of the three-dimensional integral and the surface integral; the 
sign of the second term also changes since F°° = — F". This shows that the 
Hamiltonian is not dependent on the velocities as desired. The final term, in 
which Ua is a vector normal to the surface enclosing the region being integrated, 
can be ignored if we assume that the strength of the electric field falls off at 
least as fast as 1/r^. 

Our next step is to examine the constraints. We have already found one 
constraint from the antisymmetry of the field tensor, namely that 

F°(x) fs 0. (A.45) 
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We note that this is a constraint on i5°(x) at every point in the three-space under 
consideration. To find the secondary constraint associated with this primary 
constraint, we examine the time evolution of £'°(x): 

£;"(x) w 

[E°,H] w 
dE°{:K.) dH 9£;"(x) dH 
dAo{x) dEO{x) ^ dEO{x) aAo(x) ^ 

^ ^ J Ao(x)5„i?° (x)d3x) = 



aAo(x 



The second term is zero, since the momenta are totally independent of the 
coordinates. Since the integral of the first term must be zero for a variation 
of the field at a ny po int in space x, it follows that the secondary constraint 



associated with ( A. 45 ) is 

daE^ix.) w 0. (A.46) 

The reader may note that this is Gauss' law in the absence of charge; this makes 
sense, since we have been examining the source-free versions of the Maxwell 
equations. If we examine the time evolution of this new constraint, we find that 
the equation reduces identically to zero; hence, wc have one primary constraint 
and one secondary constraint for this Hamiltonian. It can be shown that these 
constraints are both first-class. 

The total Hamiltonian of the electromagnetic field is then given by 

ffT = y" ^F-'Fab + ^E-Ea d^x - I AodaE- d^x + I v{x)E° d^x, (A.47) 

where w(x) is an arbitrary scalar function over all of the region of integration. 
To obtain the extended Hamiltonian, we add in the secondary constraints: 

He^Ht + j u{x)daEa d^x (A.48) 

where u(x) is another arbitrary scalar function. We note the similarity of this 
term to the third term of the total Hamiltonian; this suggests that Aq does not 
have any physical significance, since it can have an arbitrary coefficient added 
to it. If wc wish to discard all variables without physical significance, we can 
also drop the final term of the total Hamiltonian (involving the term E'^ — 0), 
to obtain a new Hamiltonian 

H = I ^F-'Fab + ^E''Ea+u{x)daEad'x. (A.49) 
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This Hamiltonian has the advantage of containing no reference to the unphysical 
variables Aq and while still retaining some gauge freedom (in terms of the 
function m(x).) 

Finally, we note that this expression is not as unfamiliar as it might seem. 
The first term of (A.49) can be written as 

lF'^''F,,^^iBl + Bl + Bl) (A.50) 

(where the extra factor of two comes from the fact that we have summed both 
the term F°'^Fab and F^°-Fba-) Similarly, the second term can be written as 

]^E^Ea = ]^{El + El+El) (A.51) 

Noting that daE'^ = in free space, the last term goes away, and we are left 
with the familiar expression 

= i y + d^x, (A.52) 
which is quite simply the energy stored in an electromagnetic field. 
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B Simple Spin Networks and their Values 



With the definitions and identities in Section we can create a veritable panoply 
of spin networks. More importantly, we can assign numerical values to any closed 
spin network (since if a spin network is closed, we are effectively summing over 
all the indices of its component matrices.) Some important strand identities and 
closed spin networks, along with their values, are given in we cite these 
results here: 



• The n-loop, denoted A„: 



G 



+1) 



(B.l) 



The "theta", denoted 6{a,b,c): 



a+b+c _^ iy.(^ a+b—c | ^ b+c—a | ^ c+a-b -^ | 



(B.2) 



Often, we will encounter a theta-diagram with two identical strands, and 
the third strand equal to 2. In this case, the value of the theta diagram 
becomes 



(m,m,2) = (-1) 



„+i(m + 2)(m + l) 



The "Tet" , denoted Tet 



a b e 
c d f 



2m 

, whose value is given by 



(B.3) 



Tet 



a b e 
c d f 



\z 



(-1)^(5 + 1)! 



where 



N = 



™<.<Mn.(^--.)!n,(6,-.)! 



(B.4) 



a\b\dd\e\f\ 

ai = ^{a + d + e) bi = ^{b + d + e + f) 
a2 = ^{b + c + e) 62 = 5(a + c + e + /) 

03 = 5(0 + + /) b3 = ^{a + b + c + d) 

04 = i(c + rf + /) 

m = max{ai} M = min{bj} 
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The "bubble" diagram, which is proportional to a single strand: 



_ e{a,b,c) I 

— Oaa' T a 

I 



(B.5) 



The lambda-move: 



where 



= X 



I 

ab c 



A 



ab 



(-1) 



(a(a+l)+b(&+l)-c(c+l))/2 



(B.6) 



(B.7) 



The recoupling theorem, which is perhaps the most important of these 
results: 



/ 



\ 



a b i 
c d i' 



(Bi 



where the symbol on the right-hand side is the Kauffman-Lins 6-j symbol. 
Its value is given by 



a b e 
c d f 



Tet 



a b e 
c d f 



A, 



0{a,d,e)e{b,cj) 



(B.9) 



The Kauffman-Lins 6-j symbols are related to the more commonly known 
Wigner 6-j symbols; for the details of this relation, see Appendix 

Given a 2-edge grasping one edge of a trivalent vertex, it is possible to 
"slide" it onto the other two edges, using the relation: 



r 



(B.IO) 



A network with three edges branching from the vertices of a triangle can 
be written as a multiple of a simple trivalent vertex: 



V 



Tet 


a b c 




i j k 


e 


a, b, c) 



Y 



(B.ll) 
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C Relations Between Wigner 6-j Symbols and 
KaufFman-Lins 6-j Symbols 

The Wigner 6-j symbols, used in transformations between different couplings 
of three angular momenta, are very closely related to the Kauffman-Lins 6-j 
symbols used in the recoupling theorem. Since the former of these symbols have 
been studied more thoroughly than the latter, it would be to our advantage 
to find a way to translate between them. However, there seems to be some 
confusion in the literature as to the proper translation. It can be shown, using 
the explicit formulas for the Wigner 6-j symbol in Varshalovich and the 
Kauffman-Lins 6-j symbol in Kauffman and Lins |p7| , that 

A/2 B/2 F/2 
C/2 D/2 E/2 



A B E\ ^eiB,C,E)9iA,D,E) ^_^^^e+f)/2 (^.1) 



C D F j^^ AE^/9{A,B,F)e{C,D,F) 

where the subscript KL indicates the Kauffman-Lins symbol and the W sub- 
script indicates the Wigner symbol. Note the reversal of the arguments in the 
last column. 

We will also have occasion to use the following identity, relating the Tet and 
the Wigner 6-j symbol: 



. , , , ^ Tet 
A/2 B/2 F/2 



ABE 
CDF 



C/2 D/2 E/2 \^ D, E)B{B, C, E)e{A, B, F)9{C, D, F) 



(C.2) 



In general, we will omit the KL subscript to indicate a Kauffman-Lins sym- 
bol; however, a Wigner 6-j symbol will, in this work, always have the W sub- 
script present. 
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D Diffeomorphism Invariance &: the Angle Spec- 
trum 



In our discussion of the definition and properties of spin networks, we scrupu- 
lously avoided any mention of what the exact nature of these networks might be. 
Specifically, we did not mention whether these networks are merely bookkeeping 
tools, embedded in an underlying spatial manifold and drawing on the proper- 
ties of that manifold to give values to the operators; or whether this n(rt,work 
has a deeper, more intrinsic structure. While our previous discussion proceeded 
without making this distinction, it is a natural question one might have about 
the theory; one might also wonder whether the properties of the angle operator 
presented here might shed some light on this subject. 

One of the fundamental properties of general relativity is diffeomorphism 
invariance; in other words, the equations of general relativity are invariant 
under a certain class of smooth transformations known as diffeomorphisms. To 
get a better idea of what we mean by "smooth" , suppose we have an infinitely 
stretchable, unpuncturable balloon. If we were to squash it or stretch it, we 
would have a smooth transformation of the sheet. If, however, we were to pinch 
a point on the balloon and pull it up, it would not be a smooth transformation 
at this point of the sheet. (It would, however, be a smooth transformation at 
every other point of the balloon.) 

An important property of these diffeomorphisms is that they are locally lin- 
ear; in other words, if we look at a neighbourhood around any point of the 
aforementioned balloon, we can approximate the transformation as a linear 
transformation. Moreover, this approximation can be made to arbitrary ac- 
curacy by selecting a sufficiently small neighbourhood. 

Suppose, then, that our balloon has a network drawn on it. If we examine 
a vertex of this network as we smoothly squash and stretch the balloon, we 
will see that the angles between the edges meeting at this vertex will change. 
It is evident, then, that diffeomorphisms can change the angle between two 
strands if they are embedded in an underlying manifold. However, our angle 
operator does not depend on any properties of the geometry of the strands, 
merely on their labels. This would seem to suggest, then, that the angle operator 
is diffeomorphism invariant. 

There is another, more important conclusion to be drawn from the nature 
of the angle operator. The angle spectrum for any given vertex only includes 
a finite number of angles; even if we consider the set of all vertices with a 
given valence (a countable set), we would still only have a countable number of 
possible angles in this larger set. This is at odds with the classical continuum 
model of spacetime, which predicts that within the set of all n-vertices (where 
we consider a set of n rays originating at a point in space to be an n- vertex), 
there are an uncountable number of angles. 

We could try to explain this uncountability away through diffeomorphism 
invariance; it is certainly conceivable that through adequate stretching & dis- 
torting of the surrounding space, we could turn any vertex into any other vertex; 
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at the very least, we might be able to define a countable set of equivalence classes 
of vertices, where two vertices are defined as equivalent if some diffeomorphism 
changes one vertex into the other. 

Unfortunately, this is not the case. Grot and Rovelli have shown that 
while any intersection of n lines in three-dimensional space can be deformed 
into any other intersection through arbitrary transformations of the underlying 
space, linear transformations cannot do so for an arbitrary intersection with 
n > 5. Instead, for n > 5, they show that if we define as equivalent two inter- 
sections which can be obtained from each other by a linear transformation, then 
there exist an uncountable number of equivalence classes of intersections. These 
equivalence classes can therefore be parameterized by at least one continuous 
parameter; Grot and Rovelli show that the dimension of this moduli space (i.e. 
the number of continuous parameters required) for vertices of valence n is given 

by 

5 1 

d(n) > (2n ~ b)ra ~ -rn^ - -m^ , (D.l) 

where 

V4871 - 23 - 7 



m{n) 



6 



(D.2) 



This formula yields d{n) > for n = 2,3,4, but d(5) > 2, d(6) > 4, and 
increases without bound as n — > oo. Grot and Rovelli also explicitly show (as 
an example) that d(5) is exactly equal to 2, and that the equivalence classes can 
be parameterized by the numbers 

,1,2 /1/3 

Ai = ^ and A2 = ^ (D.3) 

where we have labelled the edges {1, 2, 3, 4, 5}, and is the angle between edges 
i and j. 

This result implies that even when we take diffeomorphism invariance into 
account, there is still an uncountable spectrum of possible "angles" for vertices 
with valence greater than or equal to 5. Of course, these parameters are only 
functions of what one might normally consider to be the angles of the vertex, as 
in ( |D.3| ); nevertheless, classical vertices under diffeomorphism still possess more 
degrees of freedom in possible angle values than the quantized angle operator 
would seem to allow. 

What can we conclude from these facts? For vertices with n> 5, it should 
not surprise us that we have been reduced from an uncountably infinite number 
of possible measurements to a countably infinite number. After all, such a 
reduction occurs in many simple quantum mechanical systems, most notably in 
the measurements of angular momentum — in the classical picture, the angular 
momentum of a particle can take on any real value, while in the quantum 
picture, it must be an integer multiple of h. The existence of these classical 
parameters may also imply, however, that the angle operator is "not all there 
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is" ; if these parameters do persist in the quantum theory, then some (if not all) 
angle eigenvalues are highly degenerate. 

The n = 3 case, however, is much more surprising. In this case, the angle 
operator can still produce a countably infinite number of measurements (from 
the set of all trivalent vertices); these measurements are, as we showed above, 
diffeomorphism-invariant. However, under classical diffeomorphisms, a set of 
three vectors which span an m-dimensional space can be mapped to any other 
three vectors which also span an m-dimensional subspace. If to = 3, we have 
what Grot and Rovelli call the "non-degenerate" case; in this case, any two 
vertices arc diffeomorphism equivalent. Grot and Rovelli also cite the result of 
Arnold that this holds for the m = 2 case as well. Finally, the m = 1 case is 
trivial: the only possible classical angle measurements are 0° and 180°. Hence, 
there is only a finite amount of diffeomorphism-invariant information that can be 
extracted from a classical trivalent vertex! This is a drastic difference from the 
normal state of affairs in quantum mechanics; while quantizing a given theory 
normally restricts the amount of knowledge one can have about a system, the 
angle operator actually seems to be increasing it. This suggests that there is 
something fundamentally different going on at these low-valence vertices; the 
theory predicts that our naive continuum model of space breaks down at these 
tiny length scales. 
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